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Despite  the  substantial  theoretical  studies  of  thermal  contact  conductance  in  the  past 
decades,  the  application  of  statistical  mechanics  in  this  field  has  never  been  considered. 
The  Boltzmann  statistical  model  is  attempted  in  this  study  to  predict  the  distribution  of  true 
contact  spots  between  two  metal  surfaces  on  which  asperities  are  distributed 
homogeneously  and  isotropically.  Further,  the  relationship  between  the  thermal  contact 
conductance  and  the  relative  interfacial  pressure  is  derived.  Different  assumptions  of 
surface  asperity  deformation  mechanisms  lead  to  the  development  of  four  theoretical 
models.  Among  them,  the  most  effective  model  is  the  plastic  cone  model,  in  which 
surface  roughness  is  resembled  by  conical  asperities  and  the  load-carrying  asperities  are 
deformed  plastically.  The  numerical  results  are  compared  to  a wide  range  of  published 
empirical  data  and  a good  order-of-magnitude  agreement  is  found. 
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In  this  study,  the  conventional  assumption  of  Gaussian  summit  heights  is  shown  to 
be  unnecessary.  The  probability  density  function  of  this  kind  is  determined  naturally  by 
means  of  the  Boltzmann  model. 

Experiments  were  also  conducted  in  order  to  verify  the  accuracy  of  the  predictions 
by  the  Boltzmann  model.  Electrical  contact  conductance,  simulating  the  thermal  contact 
conductance  effect  in  a vacuum  environment,  was  measured.  The  tested  surfaces  were 
made  of  Stainless  Steel  304  and  finished  by  sandblasting.  The  experimental  results 
indicate  the  presence  of  electrically  insulating  surface  films.  The  correction  factors 
applied  to  the  predictions  are  of  the  magnitude  reported  in  the  previous  studies. 
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CHAPTER  1 
INTRODUCTION 


When  two  solid  bodies  at  different  temperatures  are  brought  into  mechanical  contact, 
a relatively  high  thermal  resistance  impedes  the  heat  flow  across  the  interface.  This 
resistance,  commonly  known  as  thermal  contact  resistance,  arises  because  the  heat-flow 
path  is  lengthened  as  a result  of  the  actual  contact  area  being  only  a small  fraction  of  the 
apparent  contact  area.  Sometimes,  the  thermal  contact  resistance  is  expressed  in  its 
reciprocal  form,  called  thermal  contact  conductance. 

Engineers  often  encounter  contact  conductance  in  a wide  range  of  thermal  systems, 
including  advanced  storage  electric  water  heaters  [1],  nuclear  reactors  [2],  high-power 
density  space  vehicles  [3],  among  others.  Recently,  thermal  contact  conductance  has  been 
found  to  play  an  important  role  in  small-scale  heat  removal  systems  such  as 
microelectronics  [4]  and  in  heat  transfer  between  superconductor  films  and  substrates  [5] . 
Thus,  the  contact  conductance  has  received  renewed  interest  and  attention. 

Due  to  the  considerable  effect  of  the  contact  conductance,  many  theoretical  and 
experimental  investigations  have  been  performed  in  the  last  five  decades.  Significant 
advances  have  been  made  in  the  development  of  theoretical  models  to  predict  the  contact 
conductance.  In  the  experimental  arena,  efforts  have  been  made  in  the  use  of  contact 
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conductance  for  thermal  rectification  [6,7]  and  in  the  use  of  metallic  coatings  to  enhance 
contact  conductance  [8,9], 

Among  all  the  existing  theoretical  models,  the  use  of  statistical  mechanics  has  never 
been  attempted  even  though  it  is  a potentially  promising  approach.  The  physical 
foundation  is  based  on  the  surface  asperities  comparable  to  molecules  in  the  traditional 
sense  of  microscopic  viewpoint.  The  focus  of  this  investigation  is  thus  directed  toward 
the  application  of  statistical  mechanics  to  the  prediction  of  contact  conductance.  Four 
different  statistical  mechanics-based  models  are  developed.  In  each  case,  assumptions  and 
detailed  derivation  are  given.  Numerical  results  show  that  the  predicted  conductance 
agrees  reasonably  well  with  a wide  range  of  published  empirical  data.  It  is  thus  concluded 
that  statistical  mechanics  is  a feasible  approach. 

The  electrical  contact  conductance  between  two  stainless  steel  surfaces  was  tested. 
Due  to  the  presence  of  surface  films,  the  measured  conductance  is  substantially  lower  than 
the  prediction.  The  role  of  surface  films  in  contact  conductance  is  discussed  in  this 


dissertation. 


CHAPTER  2 

MECHANISMS  OF  THERMAL  CONTACT  CONDUCTANCE 


The  thermal  contact  conductance  effect  is  illustrated  in  Fig.  2.1.  When  heat  flows 
steadily  across  the  joint  of  two  solid  bodies  in  mechanical  contact,  the  temperature  drops 
at  the  interface.  This  phenomenon  is  well-explained  by  the  fact  that  the  actual  contact  area 
is  exceedingly  small  as  compared  to  the  apparent  contact  area  due  to  the  presence  of 
roughness  and  waviness  of  the  engaging  surfaces.  Moreover,  the  interstitial  material,  such 
as  air,  is  a poor  conductor.  A large  portion  of  the  heat  flow  thus  converges  to  the  discrete 
solid-solid  contact  spots  as  illustrated  in  Fig.  2.2.  The  increase  of  the  heat-flow  path 
length  results  in  a thermal  resistance  across  the  interface.  Its  reciprocal  is  the  thermal 
contact  conductance,  defined  as 


There  are  three  possible  modes  of  heat  transfer  at  the  interface.  They  are  described  in  the 
sections  that  follow. 
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a) 


b) 


Fig.  2.1 


Illustration  of  thermal  contact  conductance:  a)  heat  flows  through  two 
bodies  in  contact  and  b)  the  corresponding  temperature  profile 
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Fig.  2.2 


Convergence  of  heat  flow  from  a semi-infinite  medium  to  discrete  contact 
spots 


6 


2.1  Conduction  through  True  Contact  Areas 

In  most  engineering  applications,  heat  conduction  through  the  actual  contact  spots  is 
the  dominant  mode  of  heat  transfer.  Parameters  directly  involved  in  the  thermal  analysis 
of  this  heat  transfer  consist  of  the  solid  thermal  conductivities  and  the  distribution  of  the 
true  contact  spots.  The  formation  of  contact  spots  between  two  rough  surfaces  pressed 
together  is  thus  the  focus  of  this  investigation  of  contact  conductance.  This  mode  of  heat 
transfer  will  be  elaborated  on  further  in  the  literature  review  in  Chapter  3 . 

2.2  Thermal  Radiation 


Heat  can  also  travel  by  other  means.  For  the  portions  of  the  surfaces  which  are  not 
in  physical  contact,  thermal  radiation  exchange  can  take  place.  This  mode  of  heat  transfer 
can  be  modeled  by  radiation  exchange  between  two  parallel  walls  and  the  relationship 
between  the  heat  flux  and  the  surface  temperatures  is  expressed  by 


// 


1/e,  + l/e2  - 1 


(2.2) 


°SB  — 


where 


£; 


Stefan-Boltzmann  constant 
emissivity  of  wall  i,  i = 1,2 
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Unless  the  temperature  difference  at  the  interface  is  large,  (T,  - T2)  > > 0,  the 
radiation  contribution  to  the  contact  conductance  is  negligible  for  most  practical 
applications.  Fenech  and  Rohsenow  show  that  the  ratio  of  radiation  to  conduction  through 
contact  spots  does  not  exceed  1 % when  the  interface  is  at  a mean  temperature  as  high  as 
600  °C  [10]. 


2.3  Convection  and  Conduction  through  Interstitial  Medium 


When  the  interstice  is  filled  with  a fluid,  such  as  air,  heat  may  be  transferred  by 
convection  and  conduction  in  the  interstitial  medium.  The  amount  of  convective  heat 
transfer  can  be  estimated  again  by  a two  lateral  wall  model.  As  Bejan  [11]  describes,  the 
existence  of  the  naturally  driven  fluid  circulation  between  the  hot  and  the  cold  walls  can 
be  tested  by  the  Rayleigh  number,  which  is  expressed  by 


Ra, 


s„K(Th-W‘ 


a v 

th  fl 


(2.3) 


where  ggra  = gravitational  acceleration 

pth  = coefficient  of  volumetric  thermal  expansion 
ath  = thermal  diffusivity 
vfl  = kinematic  viscosity 


For  the  case  where  a pair  of  vertical  planes  are  evenly  separated  as  illustrated  in  Fig.  2.3a, 
the  critical  condition  of  fluid  motion  is  RaL  > (L/D)4.  When  the  two  walls  are 
horizontally  positioned,  the  critical  condition  becomes  different  as  illustrated  in  Figs.  2.3b 
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a) 


v 

gravity 


> 


b) 


Ig 


RaL 


cell 


> 1708 


c) 


V 


gravity  . .j.' 


stably  stratified  fluid 


Fig.  2.3  Natural  convection  between  hot  and  cold  parallel  walls  in  three  different 
extreme  cases 
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and  3c.  For  this  system  heated  from  below  as  shown  in  Fig.  2.3b,  a Benard  convection 
flow  pattern  is  induced  when  RaL  > 1708.  Conversely,  for  the  system  heated  from  the 
upper  surface  as  shown  in  Fig.  2.3c,  the  gap  is  filled  with  quiescent,  thermally  stratified 
fluid  and  thus  convection  does  not  play  a role. 

For  the  interstitial  fluid  between  nominally  flat  surfaces  in  mechanical  contact,  the 
above  critical  conditions  are  seldom  satisfied  because  of  the  small  gap  thickness,  of  the 
order  of  1 pirn.  Therefore,  the  interstice  is  filled  with  quiescent  fluid.  Hence,  heat 
transfer  by  natural  convection  does  not  occur. 

The  gap  thickness  is  still  orders  of  magnitude  greater  than  the  mean  free  path  of  gas 
molecules.  For  instance,  the  mean  free  path  of  oxygen  at  standard  conditions  is  8 x 10'2 
jam  [12].  For  this  reason,  the  heat  transfer  through  the  interstitial  material  is  entirely  ruled 
by  conduction  under  Fourier's  law.  In  most  experimental  tests,  this  mode  of  heat  transfer 
is  commonly  eliminated  by  measuring  the  thermal  contact  conductance  in  a vacuum 


environment. 


CHAPTER  3 
LITERATURE  REVIEW 


3 . 1 Review  Papers 

The  bibliography  by  Hsieh  and  Davis  listed  most  of  the  early  publications  dealing 
with  the  thermal  contact  resistance  [13].  Subsequently,  a number  of  other  reviews  have 
appeared  in  the  literature  which  summarized  a diverse  range  of  related  research  performed 
more  recently.  Snaith  et  al.  [14]  presented  a paper  that  was  specifically  aimed  at  assisting 
designers  in  understanding  the  physical  factors  involved  by  an  overview  of  the  literature 
in  the  areas  of  contact  conductance.  The  mechanisms  of  heat  transfer  across  a metal-metal 
interface  and  the  analytical  prediction  methods  were  well  presented.  Other  discussions  in 
the  review  included  parametric  effects  and  empirical  correlations. 

Madhusudana  and  Fletcher  summarized  the  work  primarily  completed  in  the  1970s 
[15].  They  reviewed  theoretical  and  experimental  investigations  dealing  only  with  the  heat 
transfer  through  the  actual  contact  spots  and  the  interstitial  materials  by  conduction. 
Theoretical,  numerical,  and  experimental  analyses  for  constrictions  of  different  shapes 
were  considered  and  summarized  in  a table.  Correlations  of  experimental  data  and  some 
special  topics  including  stacks  of  lamination,  bolted  joints,  cylindrical  joints,  and  nuclear 
reactor  fuel  elements  were  also  covered. 
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A more  recent  review,  again  by  Fletcher,  surveyed  the  literature  published  since 
1980  [16].  In  addition  to  the  primary  areas  described  earlier,  contact  conductance  in 
microelectronics  and  biomedicine,  thermal  rectification  effect,  gas  conductance,  and 
periodic  and  sliding  contacts  were  also  considered. 

3.2  Theoretical  Models 


As  mentioned  before,  the  heat  transfer  across  contacting  surfaces  in  vacuum  is  mainly 
dominated  by  conduction  through  the  discrete  contact  spots.  The  total  resistance  can  thus 
be  represented  by  all  spot  resistances,  also  known  as  the  constriction  resistances,  in  a 
parallel  connection  as 
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TOT 


Y:_L 

' K , 

C,l 


(3.1) 


To  make  this  expression  useful,  each  constriction  resistance,  R^,  and  the  overall 
distribution  of  the  true  contacts  must  be  available.  Analytical  solutions  of  the  constriction 
resistance  have  been  successfully  developed  and  validated  by  a number  of  investigators. 
The  details  will  be  given  shortly. 

The  distribution  of  contact  spots,  another  important  piece  of  information  needed,  can 
be  obtained  graphically  by  superposing  two  recorded  profilometric  traces  as  shown  in  Fig. 


3.1  [10].  The  number  of  contacts  are  then  counted  and  the  size  of  each  contact  is 
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Fig.  3.1 


Estimation  of  the  number  and  sizes  of  contact  spots  by  a graphical  means 


13 


measured.  In  spite  of  the  feasibility,  the  process  is  extremely  laborious  and  time- 
consuming.  Preference  is  thus  given  to  mathematical  modeling. 

Three  successful  models  developed  in  the  past  will  be  discussed  in  this  section  to 
illustrate  the  different  features  that  have  been  considered  previously.  Some  of  these 
features,  including  the  surface  models  and  deformation  mechanisms,  are  applied  in 
conjunction  with  Boltzmann  statistics  to  devise  an  approach  to  predict  the  thermal  contact 
conductance.  All  these  previous  models  apply  to  contact  between  nominally  flat  surfaces, 
so  that  the  waviness  effect  is  ignored.  Their  results  will  be  presented  along  with  the 
predictions  obtained  by  the  Boltzmann  statistical  model  in  Chapter  7 for  comparison. 

3.2.1  Constriction  Resistance  due  to  Single  Contact 

The  closed  form  solution  to  the  electrical  resistance  due  to  a single  circular  contact 
spot  traces  back  to  Holm  [17]  who  used  it  to  predict  the  electrical  contact  resistance. 
Under  the  condition  where  heat  transfer  by  radiation,  convection,  and  conduction  through 
the  interstitial  materials  is  insignificant,  the  thermal  resistance,  heat  rate,  and  temperature 
gradient  are  directly  analogous  to  the  electrical  resistance,  current,  and  potential  gradient, 
respectively.  Therefore,  it  is  justified  to  apply  Holm's  work  to  the  study  of  thermal 
contact  conductance. 

Holm's  constriction  resistance  can  be  derived  by  using  a Laplace  equation  in  a semi- 
infinite  medium  of  constant  properties  with  the  boundary  conditions  given  in  Fig.  3.2. 
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Fig.  3.2 


Heat  flows  through  a single  constriction  in  a semi-infinite  medium 
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The  resultant  resistance  is  obtained  as 


R 

C 


1 

4a  k 


(3.2) 


When  the  resistance  is  across  two  semi-infinite  media  of  different  materials,  as  illustrated 
in  Fig.  3.3,  the  contact  surfaces  may  be  considered  isothermal.  Then,  the  total  resistance 
can  be  represented  as  a sum  as 


For  similar  materials, 


thus, 
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1 
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(3.3) 


(3.4) 


The  above  model  provides  accurate  results  when  the  actual  area  of  contact  is  only  a 
small  fraction  of  the  apparent  area.  As  this  fraction  increases,  heat  traveling  through  a 
contact  spot  from  a semi-infinite  cylindrical  medium  becomes  more  realistic  in  analyzing 
the  heat  transfer.  Figure  3.4  shows  the  schematic  of  this  cylindrical  heat  channel, 
including  the  boundary  conditions.  The  analytical  solution  to  this  system  leads  to  a 
correction  of  Eq.  3.4  by  a multiplying  factor,  namely,  the  constriction  alleviation  factor, 
Y,  as 


2 ak 


(3.5) 
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Fig.  3.3 


Heat  flows  across  two  semi-infinite  media  through  a single 
constriction 


17 


31 

dz 


Fig.  3.4 


Cylindrical  flow  channel  and  the  corresponding  boundary  conditions 
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A number  of  investigators  have  successfully  derived  analytical  and  approximate  solutions 
to  this  model.  Gibson  reduced  the  mixed  boundary  value  problem  to  a Fredholm  integral 
equation  and  then  solved  it  [18].  A similar  expression  was  obtained  by  Roess  as  presented 
by  Snaith  et  al.  [14].  Approximations  given  by  Cooper  et  al.  [19]  were  rather  simple  to 
use.  All  results  mentioned  above  are  summarized  in  Table  3.1. 

3.2.2  Hsieh  Model 


Incorporating  the  constriction  resistance  analysis  given  above,  Hsieh  performed  a 
statistical  analysis  to  predict  the  total  contact  conductance  [20] . His  model  is  constructed 
on  the  basis  of  contact  between  a pair  of  homogeneous  and  isotropic  surfaces.  The 
irregularities  of  each  rough  surface  are  modeled  by  conical  asperities  of  an  identical 
gradient  and  various  heights  resting  on  a common  flat  plane;  see  Fig.  3.5.  This  surface 
model  has  been  proven  to  be  a suitable  representation  of  the  surfaces  processed  by 
grinding,  lapping,  honing,  bead-blasting,  anodizing,  and  accurate  casting  [21,22],  In 
addition,  the  summit  height  distribution  is  assumed  Gaussian. 

From  the  mechanical  deformation  point  of  view,  touching  asperities  resemble  blunt 
indenters.  This  results  in  a plastic  summit  deformation  and  a local  elastic  substrate 
deformation  as  shown  in  Fig.  3.6.  Working  collectively,  they  serve  well  to  estimate  the 
load  distribution  over  the  apparent  area. 

By  using  the  above  assumptions  and  a given  base  plane  separation  which  is  small 
enough  to  initiate  contact,  the  probability  density  of  surface  asperities  making  contacts  is 


Table  3.1 


Summary  of  constriction  alleviation  factors  obtained  by  Gibson,  Roess,  and 
Cooper  et  al. 


Gibson  [18] 

* = 1 - 1.4091839(a/b)  + 0.338010(a/b)3  + 0.067902(a/b)5 
+ 0((a/b)6) 

Roess  [14] 

* = 1 - 1. 40925  (a/b)  + 0.29591(a/b)3  + 0.0525419(a/b)5 
+ 0((a/b)6) 

Cooper  et  al. 
[19] 

/ (1  - a/b)1 5 , all  a/b 

¥ = J 1 - 1.4(a/b)  , a/b  < 0.5 

l 1 , small  a/b 
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Fig.  3.5 


Surface  roughness  modeled  by  conical  asperities 
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Fig.  3.6 


Asperity  deformation  mechanism  used  by  the  Hsieh  model  [20] 
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derived.  Knowing  this  probability  density  function,  Hsieh  finally  evaluated  the 
"expectations"  of  the  total  thermal  contact  conductance  and  the  relative  total  load. 

3.2.3  Yovanovich  Model 

In  the  past,  a number  of  correlations  have  been  proposed.  The  Yovanovich  thermal 
contact  model  is  one  of  the  successful  examples.  This  model  was  first  presented  by 
Cooper  et  al.  [19]  and  later  refined  by  Yovanovich  [23].  Unlike  the  Hsieh  model,  it  is 
assumed  that  the  surface  heights  (not  summit  heights)  form  a Gaussian  distribution. 
Moreover,  the  mean  surface  profile  slope  m is  independent  of  the  surface  height.  The 
contact  between  two  rough  surfaces  is  interpreted  by  a contact  between  a perfectly  flat 
surface  and  an  equivalent  rough  surface  having  combined  surface  parameters  as  illustrated 
in  Fig.  3.7.  These  parameters  lead  to  a combined  RMS  roughness,  defined  as 


When  two  surfaces  are  pressing  together,  it  is  assumed  that  the  summits  of  all 


(3.6) 


and  a combined  mean  slope,  defined  as 


(3.7) 


contacting  asperities  are  deformed  plastically  while  the  entire  substrate  undergoes  an 
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Fig.  3.7 


a)  Contacting  rough  surfaces  modeled  by  b)  the  contact  between  a flat 
plane  and  an  equivalent  rough  surface 
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elastic  deformation  (absence  of  local  substrate  deformation).  By  coupling  this  deformation 
mechanism  with  the  surface  statistical  properties  and  the  constriction  resistance  equation 
(Eq.  3.4),  an  expression  giving  a relationship  among  the  contact  conductance,  the  contact 
pressure,  and  other  parameters  involved  can  be  derived.  This  expression  is  then  correlated 
by  a simple  equation  as 


h o 
km 


1.25 


0.95 


H 


(3.8) 


where  H and  k are  the  harmonic  mean  hardness  and  conductivity,  respectively.  If  the 
surfaces  in  contact  are  made  of  dissimilar  materials, 

J_ 

H 


k 

Due  to  its  accuracy  and  simplicity,  the  Yovanovich  model  is  widely  accepted  and  utilized. 

3.2.4  Modified  Greenwood-Williamson  Model 

With  the  belief  that  both  plastic  and  elastic  deformations  of  the  contact  asperities  take 
place  in  practice,  McWaid  and  Marschall  [24]  investigated  an  elastic  contact  model. 
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Known  as  the  modified  Greenwood-Williamson  model,  the  method  owes  its  origin  to  the 
work  initially  performed  by  Greenwood  and  Williamson  [25]. 

The  assumptions  made  here  are  quite  different  from  the  ones  described  previously. 
Since  the  summit  of  an  asperity  is  commonly  round,  a rough  surface  is  postulated  as  a flat 
plane  covered  with  spherical  asperities  of  an  identical  mean  curvature  as  shown  in  Fig. 
3.8.  The  fundamental  surface  parameters  are  the  zeroth,  second,  and  fourth  moments  of 
the  power  spectral  density,  denoted  by  mo0,  mo2,  and  mo4,  respectively.  Several 
topographic  parameters,  such  as  the  standard  deviation  of  the  summit  height,  the  density 
of  summits,  etc.,  can  be  expressed  in  terms  of  these  spectral  moments. 

The  estimation  of  the  degree  of  engagement  of  the  contact  is  treated  in  a similar 
fashion  as  Yovanovich.  The  contact  between  two  rough  surfaces  is  represented  by  the 
contact  between  an  ideally  smooth  plane  and  an  equivalent  rough  surface  whose  spectral 
moments  equal  to  the  sum  of  the  respective  moments  of  the  individual  surfaces. 
Additionally,  this  model  assumes  that  the  summit  heights  of  the  equivalent  rough  surface 
form  a Gaussian  distribution.  The  elastic  deformation  analysis  of  a single  contact  is 
accomplished  by  invoking  the  use  of  the  Hertzian  theory  for  the  contact  between  a sphere 
and  a plane  [26] . 

Using  the  assumptions  and  analyses  given  above,  one  can  obtain  the  expressions  of 
the  loading  pressure  P and  the  conductance  h in  terms  of  the  mean  plane  separation  d as 


(3.9) 
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Fig.  3.8 


Surface  roughness  modeled  by  spherical  asperities 
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Hence,  a relationship  between  h and  P is  established. 

McWaid  and  Marschall  assume  that  all  contacts  undergo  an  elastic  deformation.  On 
the  other  hand,  Tabor  suggests  a criterion  for  the  onset  of  plastic  deformation  of  a 
spherical  asperity  [27] . It  is  shown  that  the  first  plastic  flow  occurs  when  the  maximum 
Hertzian  pressure  is  greater  than  60%  of  the  indentation  hardness.  In  other  words,  if  the 
maximum  pressure  is  less  than  that,  the  entire  asperity  is  deformed  elastically.  Both 
mechanisms  of  deformation  of  the  spherical  asperities  will  be  employed  in  conjunction 
with  the  Boltzmann  statistical  model  in  analyzing  the  contact  conductance  in  Chapter  6. 
More  details  of  this  modeling  technique,  including  the  topographic  parameters  and  the  use 
of  the  Hertzian  theory,  will  be  discussed  then. 
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3.3  Experimental  Studies 


3.3.1  Experimental  Setup 

In  addition  to  the  theoretical  studies,  numerous  efforts  have  been  made  to  measure 
the  contact  conductance  experimentally.  Mainly  for  the  purpose  of  validating  analytical 
models,  the  experimental  measurements  are  indispensable  in  assessing  the  effects  of  the 
contacting  materials,  the  interstitial  fluids,  or  other  characteristics  of  the  surface  on  the 
contact  conductance. 

The  schematic  of  an  apparatus  setup  that  appears  repeatedly  in  the  literature  is  shown 
in  Fig.  3.9.  The  tested  samples  are  drilled  sideways  for  installation  of  thermocouples. 
After  the  processed  surfaces  are  cleaned  and  their  topographic  parameters  are  measured 
properly,  the  surfaces  are  stacked  together  between  a heat  source  and  a heat  sink.  The 
role  of  the  ball  bearings  is  to  minimize  the  effect  due  to  misalignment.  The  radiation 
shield  is  employed  to  reduce  the  radial  heat  transfer  from  the  test  specimens  by  radiation. 
Thus,  in  a vacuum  environment  heat  is  only  allowed  to  flow  axially  through  the 
specimens.  As  a uniform  heat  flux  is  introduced  by  the  heat  source,  the  temperatures  at 
the  interface  between  the  two  samples  in  contact  are  determined  by  extrapolating  the 
temperature  measurements  in  the  specimens.  The  temperature  drop,  AT,  is  then 
determined  with  the  facility  and  the  thermal  contact  conductance  can  be  obtained  by  Eq. 
2.1.  Miller  and  Fletcher  [28]  explain  each  major  component  of  this  measuring  facility  in 


more  detail. 
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Schematic  of  experimental  test  facility 
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3.3.2  Sources  of  Experimental  Data 

A large  number  of  references  of  experimental  measurements  of  thermal  contact 
conductance  can  be  found  in  the  literature  where  correlations  are  given.  Tien  [29]  listed 
twelve  different  experimental  studies  on  which  correlation  was  performed.  They  included 
tests  of  aluminum  and  stainless  steel  specimens  in  a vacuum  environment.  References  to 
tests  under  the  presence  of  various  kinds  of  interstitial  fluids,  such  as  air,  oil,  argon, 
helium,  and  paraffin,  were  provided  in  Veziroglu's  correlation  study  [30], 

A comprehensive  test  of  the  thermal  contact  conductance  was  recently  conducted  by 
Hegazy  [31].  The  surface  properties  were  characterized  in  detail  and  the  specimens  tested 
include  Nickel  200,  Stainless  Steel  304,  Zircaloy-4,  and  Zirconium-2.5  %-Nb.  The 
experiments  were  carried  out  in  vacuum,  nitrogen,  and  helium  environments. 

Surface  roughness  is  commonly  identified  by  the  measurements  of  standard  deviation 
of  the  surface  height,  o,  and  the  mean  surface  slope,  m.  However,  the  surface  slope  was 
not  commonly  reported  in  most  early  experimental  works  so  that  the  use  of  their  data  is 
rather  limited.  A correlation  performed  by  Antonetti  and  Whittle  [32]  is  worth  noting. 
They  gathered  published  topographic  data  for  65  specimens  to  establish  the  relationship 
between  the  average  roughness  and  the  root-mean-square  asperity  slope.  The  resultant 
correlation  can  also  be  expressed  in  terms  of  the  standard  deviation  of  the  surface 
roughness,  o,  in  micrometers  and  the  mean  asperity  slope,  m,  by 


m = 0.124  a 0743 


(3.11) 
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For  the  modified  Greenwood  and  Williamson  model,  in  which  spherical  asperities  are 
used,  o and  m are  not  sufficient  to  characterize  the  surface  topography.  The  three  spectral 
moments  mo0,  mo2,  and  mo4  are  the  fundamental  parameters  needed.  At  present,  the  only 
investigation  in  which  both  the  spectral  moments  and  the  thermal  contact  conductance  data 
are  reported  is  the  reference  by  McWaid  and  Marschall  [24]. 

The  numerical  results  by  the  Boltzmann  statistics-based  models  are  compared  to  the 
data  reported  in  some  of  the  references  quoted  above.  Only  those  experiments  conducted 
in  vacuum  will  be  compared  as  shown  in  Chapter  7. 


CHAPTER  4 

THEORETICAL  BASIS  FOR  THE  USE  OF  STATISTICAL  MECHANICS 


Statistical  mechanics  has  never  been  attempted  in  the  prediction  of  thermal  contact 
conductance.  In  fact,  the  problem  domain  of  the  contact  conductance  can  be  readily 
explored  by  means  of  statistical  mechanics.  The  purpose  of  this  chapter  is  to  demonstrate 
this  applicability  from  the  microscopic  point  of  view. 

In  principle,  if  an  observable  property  of  a many  particle  system  at  a state  of 
equilibrium  can  be  characterized  by  a small  number  of  other  macroscopic  properties,  this 
observable  property  can  often  be  described  by  statistical  mechanics.  This  analysis  also 
provides  a microscopic  description  of  the  system  (distribution  of  particles).  Frequently, 
statistical  mechanics  is  illustrated  by  applying  it  to  equilibrium  thermodynamic  quantities. 
An  illustrative  example  in  statistical  thermodynamics  is  further  introduced  in  this  chapter 
in  order  to  establish  the  explicit  links  between  thermal  contact  conductance  and  statistical 
mechanics. 

4.1  Microscopic  and  Macroscopic  Points  of  View  in  Contact  Conductance 

The  mechanical  and  thermal  mechanisms  for  metallic  surfaces  in  contact  in  vacuum 
have  been  discussed  in  Chapter  2.  The  macroscopic  or  observable  properties  involved  in 
this  system  are  mainly  the  loading  interfacial  pressure  and  the  thermal  contact 
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conductance.  As  one  studies  the  system  microscopically,  it  is  rather  obvious  that  the 
applied  load  at  the  interface  is  distributed  over  the  contact  spots  formed  by  the  surface 
asperities  (analogous  to  particles).  These  asperities  on  opposing  surfaces  are  first  brought 
into  contact  and  they  are  deformed  until  they  reach  a state  of  force  "equilibrium. " 

A thermal  contact  resistance,  an  observable  property,  arises  as  the  heat  is  limited  to 
flow  through  these  contact  spots.  Substantial  experimental  studies  have  shown  that,  for 
surfaces  of  given  material  properties  and  topographic  characteristics,  the  reduction  of 
contact  resistance  behaves  orderly  and  consistently  as  the  loading  pressure  increases. 
Therefore,  it  is  believed  that  the  distribution  of  the  contact  spots  is  highly  repeatable  when 
the  surfaces  are  brought  into  contact  under  a specific  loading  pressure  condition.  In  other 
words,  the  contact  conductance  can  be  considered  a result  of  "the  most  probable 
distribution"  of  contacts.  From  this  perspective,  statistical  mechanics  is  indeed  an 
appropriate  method  to  predict  the  thermal  contact  conductance. 

4.2  Boltzmann  Statistics  and  Thermodynamics 

Boltzmann  statistical  model  is  used  repeatedly  in  this  study  to  obtain  first  the  most 
probable  distribution  of  the  surface  summit  heights,  which  is,  in  turn,  used  in  the  analysis 
of  the  distribution  of  contact  spots  at  various  loading  pressures.  This  statistical  model  is 
briefly  introduced  here  from  a statistical  thermodynamics  viewpoint.  The  detailed 
discussion  can  be  found  in  the  reference  by  Sonntag  and  Van  Wylen  [33] . 
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Consider  a system  comprised  of  N particles  which  are  independently  distributed 
among  various  energy  levels.  At  the  energy  level  up  there  is  a corresponding  statistical 
weight  gj,  which  is  the  number  of  quantum  states  having  that  energy.  For  Boltzmann 
statistics,  it  is  assumed  that  all  these  particles  are  distinguishable  and  there  is  no  limit  on 
the  number  of  particles  per  quantum  state.  Then,  the  number  of  ways  that  a given 
distribution  can  arise  is  given  by  the  formula 


where  Nj  is  the  number  of  particles  having  energy  ur 

The  quantity  w given  above  is  commonly  known  as  the  thermodynamic  probability. 
The  product  n can  be  removed  by  taking  logarithm  of  both  sides  of  the  equation,  giving 


The  two  constraints  which  govern  the  activity  of  the  particles  in  a closed  system  under  an 
equilibrium  state  in  thermodynamics  are  1)  the  total  number  of  particles  must  be 
conserved,  i.e., 


(4.1) 


r'N.  = N 


(4.2) 
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and  2)  the  internal  energy  must  be  conserved,  i.e.. 


N.u.  = U 
J ' ' 


(4.3) 


To  find  the  most  probable  distribution  of  particles  NJ5  which  is  the  distribution  that  gives 
the  highest  thermodynamic  probability,  the  following  mathematical  procedures  are 
performed. 

First,  differentiate  Eq.  4.1  and  set  it  to  zero.  It  follows  that 

rf(lnw)  = ^(ln g.dN.  - In N.dN.  - dN.)  = 0 (4.4) 


Furthermore,  since  N and  U in  Eqs.  4.2  and  3 are  constants,  the  derivatives  of  these 
equations  give 


? dN.  = dN  = 0 
J 1 


(4.5) 


and 


? u.dN.  = dU  = 0 
J 1 1 


Equation  4.4  can  thus  be  simplified  by  using  Eq.  4.5  as 


(4.6) 


<i(lnw)  = 


E 

j 


In 


{ \ 

A 

\NU 


dN.  = 0 

j 


(4.7) 
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At  this  point,  the  expression  giving  the  maximum  of  w can  be  found  by  the  method  of 
Lagrange  multipliers.  Equations  4.5  and  6 are  multiplied  respectively  by  arbitrary  values 
a and  p and  added  together.  Then,  subtracting  Eq.  4.7  from  the  sum  gives 


E 

j 


( 

\ 


N. 

In—  + a 


\ 


dN.  = 0 

J 


(4.8) 


The  only  way  that  Eq.  4.8  can  be  satisfied  is  for  each  of  the  coefficients  to  be  identically 
zero,  or 

N. 

In—  + a + Bm.  = 0 

Alternatively,  it  can  be  written  in  the  form  commonly  known  as  the  Boltzmann  distribution 
law 

N.  = g.e  ae  (4.9) 

/ °i 


Finally,  the  above  equation  can  be  further  simplified  by  eliminating  a by  Eq.  4.2 
(Constraint  1).  Hence,  the  most  probable  distribution  is  obtained  as 


N. 

j 


Ngse "p“' 


•j 


(4.10) 
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When  the  discrete  variables  u,  and  g-}  are  so  closely-spaced  as  to  form  a continuum, 
the  summation  is  changed  to  an  integration.  This  is  a transformation  means  of  the  Euler- 
Maclaurin  summation  formula.  Consequently,  Eq.  4.10  is  changed  to  the  following 


In  practice,  the  computation  of  the  integral  can  be  accomplished  either  analytically  or 
numerically. 

4.3  Comparison  between  Contact  Conductance  and  Statistical  Thermodynamics 

Attention  is  now  turned  to  the  comparison  of  the  features  associated  with  the  contact 
conductance  and  the  counterpart  in  statistical  thermodynamics  given  above  to  demonstrate 
their  similarities.  They  are  identified  and  listed  side-by-side  in  Table  4.1.  This  table 
should  provide  the  rationale  of  determining  the  contact  conductance  by  statistical 
mechanics. 

Undoubtedly,  due  to  the  similarities  provided,  the  mathematics  given  in  this  chapter 
is  to  be  performed  as  the  Boltzmann  statistical  model  is  applied  to  the  prediction  of  the 
contact  conductance.  Yet,  an  appropriate  statistical  weight  function  has  not  been  defined. 
Chapter  6 will  give  the  details  of  the  derivation  of  four  statistical  weights  and  their 


dN 

U 


Ndu 


(4.11) 


o 


formulated  distributions  of  contacts. 
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Table  4. 1 Similarities  between  the  domain  of  statistical  thermodynamics  and  that  of 
thermal  contact  conductance 


Statistical  Thermodynamics 

Thermal  Contact  Conductance 

1 

Gas  particles  are  independent  of  each 
other 

Asperities  are  independent  of  each 
other 

2 

Number  of  particles  is  conserved  in 
a closed  system, 

Number  of  asperities  present  on  each 
surface  is  conserved. 

n 

2*  Ni  = N 
J 3 

3 

Total  energy  is  conserved, 

S N.  u , = U 
3 3 3 

Sum  of  loads  carried  by  asperities 
equals  the  applied  force, 

E ^ F.  = P 
3 3 3 

4 

Most  probable  distribution  of 
particles  is  found  by  the  Boltzmann 
statistical  model;  other  system 
characteristics  can  be  determined 
accordingly 

Most  probable  distribution  of  true 
contact  areas  is  found  by  the 
Boltzmann  statistical  model;  thermal 
contact  conductance  can  be 
determined  accordingly 

CHAPTER  5 

MOST  PROBABLE  DISTRIBUTION  OF  SUMMIT  HEIGHTS 
BY  BOLTZMANN  STATISTICAL  MODEL 


Before  studying  the  interaction  between  two  contacting  surfaces,  one  should  identify 
the  statistical  features  of  each  surface.  Indeed,  statistical  mechanics  can  play  a vital  role 
in  determining  the  distribution  of  the  summit  heights,  which  is  commonly  assumed 
Gaussian. 


5.1  Surface  Measurement  as  Constraints 

The  topography  of  a surface  is  often  characterized  by  the  use  of  a profilometer.  As 
the  stylus  of  a profilometer  travels  along  a straight  line  on  a surface,  the  electro- 
mechanical transducer  produces  electrical  signals  proportional  to  the  variation  of  the 
surface  heights.  The  actual  surface  profile  can  also  be  drawn  by  a chart  recorder.  For 
some  surface  preparation  techniques,  such  as  sandblasting,  the  profilometer  tends  to 
measure  a consistent  value  of  standard  deviation. 

Still,  the  above  measurements  cannot  furnish  the  detailed  statistics  of  the  summit 
heights  because  the  stylus  trace  can  pass  over  an  asperity  shoulder  rather  than  its  summit. 
The  variation  of  the  summit  heights  is  rarely  measured  due  to  the  difficulty  in  setting  up 
the  instruments.  Nevertheless,  Williamson  [34,35]  implemented  an  unconventional 
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scheme  to  examine  several  surface  texture  parameters,  including  the  summit  heights.  This 
technique  is  simple  to  understand  but  hard  to  perform.  In  principle,  this  method  creates 
a contour  map  of  a surface  by  displaying  many  closely-spaced  parallel  profiles  in  a two- 
dimensional  array.  Then,  the  summit  of  each  asperity  can  be  located.  The  major 
difficulty  in  this  method  is  in  maintaining  a unique  reference  height  of  all  profiles.  This 
problem  is  solved  by  an  auxiliary  stylus-lifter  built  by  Williamson. 

One  of  the  examples  studied  by  Williamson  is  a bead-blasted  aluminum  surface, 
whose  texture  is  highly  homogeneous  and  isotropic.  Both  the  surface  heights  and  the 
summit  heights  on  this  surface  are  reported.  The  data  show  that  there  is  a constant  shift 
S on  their  distribution  curves  as  shown  in  Fig.  5.1.  This  result  gives  the  indication  that 
the  summit  heights  and  surface  heights  should  have  the  same  standard  deviation, 

o . = O . = 0 

asperity  surface  measured 

and  their  mean  values,  fj.,  differ  by  the  value  of  S, 

^ asperity  ^ surface  ^ 

The  discussion  given  above  suggests  that  a bead-blasted  surface  possesses  the 
following  characteristics: 

1)  Density  of  asperities, 

= N 

] 1 


(5.1) 


F reguency 
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Fig.  5.1 


Distributions  of  surface  heights  and  summit  heights  of  a bead-blasted 
surface 
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2)  Mean  summit  height, 


? ■ we 


(5.2) 


3)  Standard  deviation  of  summit  heights, 


E 

j 


N.(x.  - M)2 


No2 


(5.3) 


where 


N = density  of  asperities  (number  per  unit  area) 
Nj  = density  of  asperities  having  summit  height  Xj 
Xj  = particular  summit  height 


These  relationships  can  serve  as  the  constraints  imposed  on  the  surface  topography. 


5.2  Derivation  by  Boltzmann  Statistical  Model 


The  Boltzmann  statistical  model  is  used  to  find  the  most  probable  distribution  of 
summit  heights.  The  mathematical  procedures  given  in  Section  4.2  are  used.  The 
derivatives  of  the  constraints  are: 


1) 


? dN  = 0 
J 1 


(5.4) 


2) 


? x.dN.  = 0 

J J 1 


E 

j 


(Xj  - n)2  dN. 


0 


(5.5) 


3) 


(5.6) 
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Equation  5.6  can  be  expanded  into 

XdN.  = 0 
J ' 

The  second  and  the  third  terms  have  been  shown  to  be  zero  by  Eqs.  5.5  and  4, 
respectively.  Thus,  the  derivative  of  the  third  constraint  becomes 


^ xf  dN.  - 2 n ? x dN.  + u2 
J j j J 1 1 * 


? xfdN.  = 0 
J 1 1 


(5.7) 


Applying  multipliers  a,  p,  and  y to  Eqs.  5.4,  5,  and  7,  respectively,  and  following 
the  mathematics  given  in  Eqs.  4.7,  8,  and  9,  one  can  obtain  the  Boltzmann  distribution  as 


N. 

j 


*j' 


M, 


y*i 


(5.8) 


Substituting  Eq.  5.8  into  Eq.  5.1  yields 


N 


E 

j 


p x,  -yx. 

gj  e ‘e  > 


Thus,  the  multiplier  a can  be  eliminated  in  Eq.  5.8  and  the  distribution  becomes 


N.  = 

j 


Ng.e*x‘eyx‘ 


gie 


. P _ -yx. 


(5.9) 
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Substituting  this  Nj  into  constraints  2 and  3 (Eqs.  5.2  and  3)  yields 


- p x,  yx, 

gj  e ' e 1 x. 


Px.  -yx, 

gj  e ‘e  1 


(5.10) 


and 


^ 8je'fXte'1Xi(Xj  - n)2 

J- - . o1  (5.11) 

Further,  applying  Euler-Maclaurin  summation  formula  to  convert  the  above  discrete  forms 
to  continuous  forms  yields 


dN 

X 

~Ndx 


J gxe  Px e'yx2dx 


(5.12) 


j gxe  *x e yx2xdx 

= /x  (5.13) 

J gx  e p*  e yx*dx 
0 

and 


J gx  e Pjc  e yx\x  - n)2  dx 


J gx  e Px  e yx" dx 
0 


(5.14) 
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As  the  summit  height  x is  statistically  weighted  evenly,  gx  is  a constant.  Then,  Eqs. 
5.12,  13,  and  14  become 


dNx  e e _Y*2 
Ndx 

J e f>xe  yxldx 


(5.15) 


Je  *x  e yx\dx 

= m (5.16) 

J e e ~y%1  dx 


and 


j' e *x  e yx\x  - n)2  dx 

- — = o2  (5.17) 

J e-»x  e yxldx 


Here,  the  two  unknown  multipliers,  p and  y,  are  obtained  by  solving  Eqs.  5.16  and  17. 
First,  the  integral  in  the  denominator  is  performed  exactly  as 


/ 


e pjt  e ~yx  dx 


I — JL 

V'TC  p 4y 

2^ 


lim 


erf 


\fyx 


2\/y, 


erf 


2\fy 


(5.18) 
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At  this  point,  two  hypotheses  are  made  to  simplify  the  derivation  and  their  validations 
are  provided  later.  Hypothesis  1 is  given  as 


erf 


2/y 


* -1 


(5.19) 


Thus,  Eq.  5.18  becomes 


/ 


e~Px  e y*2dx  = 


\ 


TC 

— e 


4y 


(5.20) 


Using  integration  by  parts  and  Eq.  5.20,  one  can  solve  the  numerators  in  Eqs.  5.16  and 
17  as 


' 

r2> 

f e Px  e y%l xdx  = — — 

-1  ♦ pA 

p 

— e 4y 

J 2 v 

0 1 

l \ 

y y 

and 


J~e  Pxe  yx\x  n)2dx 


— —-2  ix 
2Y 


W 1N 


(5.22) 


V 2 Y/ 


2 Y 


-2/1 


- * - 
2 Y / V 2y 


J\ 


\ 


TC 


4y 


y y 


Substituting  Eqs.  5.20,  21,  22  into  Eqs.  5.16  and  17  yields 


1 

2y 


A 


7T 

— e 


2y 


0 


4y 


(5.23) 
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and 


-&-  + 2 n 

2y  . 


2 Y 


— <?  4* 


y-  + 2M 

2y 


_P 
2 y y 


— -M2 

2 Y 


=0  (5.24) 


Hypothesis  2 is  stated  that  the  polynomial  terms  are  dominant  so  that  the  exponential 
terms  are  neglected  in  Eqs.  5.23  and  24.  Thus,  the  two  equations  become 


♦ n = 0 

2 Y 


(5.25) 


and 


2 Y 


2/T 


/ J_N 

\ 2Y/ 


— — /i2|  - o2  = 0 
2 Y 


(5.26) 


The  above  two  equations  are  solved  for  the  values  of  p and  y as 


and 


(5.27) 


Y = 


2o2 


(5.28) 


The  two  hypotheses  made  are  now  to  be  validated.  The  values  of  p and  y are 
substituted  into  the  left  hand  side  of  Eq.  5.19.  This  results  in 


erf 


1 P ^ 

V Vy/ 


erf 


l \ 


V 


s/2o 
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Surface  measurement  has  shown  /x  is  approximately  4 o [20],  thus 


erf 


/ 

V 


2Vy/ 


erfir  2.9294) 

- 0.99993 

- 1 


Hence,  Hypothesis  1 is  validated. 

The  validation  of  Hypothesis  2 is  accomplished  by  applying  order-of-magnitude 
analyses  on  Eqs.  5.23  and  24.  The  results  are  as  follows: 

1)  In  Eq.  5.23,  the  analysis  shows 


1 

l 

p 

2 Y 

p2 

l 2 Y 

—e  4y 
Y 

o 

\fTiie  2°2 


+ 0*  - /*)  =o 


order  of 


( \ 

a 


e2°  I 
\ e 


< < order  of  /x 


2)  In  Eq.  5.24,  the  analysis  shows 


-&-+2 n 

2 y 1 


2 Y 

P2 

\ 

— e 4y 
Y 

(-#* 

+ 2fi)o2 

s/2 

/ 

2y 


_P_' 

2 Y / 


2 Y 


- o 


= 0 


[(-p  + 2fi)(-fx)  + (o2  + jx2)  - oz]  =0 


20' 


order  of 


HQ 


a 2o 
\ e 


< < order  of  /x2 


Hence,  Hypothesis  2 is  validated. 


49 


Finally,  by  substituting  the  values  of  p,  y,  and  Eq.  5.20  into  Eq.  5.15,  the  most 
probable  distribution  of  summit  heights  is  expressed  by 


dN 

X 

~Ndx 


fix) 


(5.29) 


It  is  found  that  the  above  probability  density  function  for  the  Boltzmann  distribution  is 
identical  to  the  one  for  a Gaussian  distribution.  Hence,  statistical  mechanics  has  provided 
an  analytical  foundation  to  support  the  previously  made  assumption  of  Gaussian 
distribution  of  summit  heights  by  Hsieh  [20]  and  McWaid  and  Marschall  [24] . 


CHAPTER  6 

PREDICTION  OF  THERMAL  CONTACT  CONDUCTANCE  BY 
BOLTZMANN  STATISTICAL  MODEL 


In  the  preceding  chapters,  the  important  components  of  the  thermal  contact 
conductance  and  the  Boltzmann  statistical  model  have  been  reviewed.  Using  these  as 
background,  four  predicting  models  will  be  developed. 

As  reviewed  previously,  alternative  surface  models  (conical  and  spherical  asperities) 
and  deformation  mechanisms  (elastic  and  plastic)  have  been  assumed  by  different 
investigators  for  the  prediction  of  thermal  contact  conductance  with  reasonably  accurate 
results.  Thus,  for  a thorough  exploration,  the  Boltzmann  statistical  model  is  applied  in 
conjunction  with  the  combinations  of  the  above  assumptions.  Four  different  Boltzmann 
statistics-based  models  are  then  developed.  This  investigation  attempts  to  identify  the  best 
surface  model  and  asperity  deformation  mechanism  in  the  statistical  mechanics  approach. 

For  simplicity,  all  models  are  set  up  on  the  basis  of  contacts  of  similar  material  and 
surface  roughness.  It  is  expected  that  the  analytical  models  can  be  extended  for  the 
prediction  of  the  conductance  of  contacts  between  dissimilar  materials. 
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6.1  Development  of  Core  Equations 

Knowing  the  statistical  distribution  of  summit  heights  on  a surface,  derived  in 
Chapter  5 , one  can  predict  both  the  microscopic  and  macroscopic  properties  of  surfaces 
in  contact  by  the  Boltzmann  model.  The  formation  of  true  contact  areas  are  modeled  as 
demonstrated  in  Fig.  6.1,  where  cones  are  selected  for  the  shape  of  the  asperities  for 
illustration  purpose.  Figure  6.  la  shows  that  the  contact  spots  often  occur  near  the  summits 
as  Williamson  [34]  found  in  his  study.  The  first  step  in  modeling  is  to  represent  each 
rough  surface  by  the  ideally  shaped  asperities  as  illustrated  in  Fig.  6.1b.  As  expected, 
perfect  alignment  between  a pair  of  touching  asperities  rarely  occurs.  Nevertheless,  it  is 
believed  that  the  degree  of  misalignments  is  statistically  constant  since  the  macroscopic 
properties,  such  as  contact  conductance,  surface  friction,  etc.,  are  consistently 
reproducible.  It  is  thus  justified  to  simplify  the  analysis  by  assuming  that  all  touching 
asperities  are  perfectly  aligned  as  illustrated  in  Fig.  6.1c.  The  theoretical  prediction  of  the 
contact  conductance  can  be  corrected,  if  necessary,  by  a multiplying  factor,  which  is  close 
to  unity  [20] . 

Having  the  contact  modeled  as  described,  one  is  ready  to  apply  the  Boltzmann  model 
to  determine  the  distribution  of  contact  spots.  Intuitively,  the  size  of  a contact  can  be 
quantified  by  its  radius.  Alternatively,  it  can  be  represented  implicitly  by  the  summit 
penetration  depth  6 (see  Fig.  6. 1).  The  latter  is  preferred  because  the  Boltzmann  statistical 
model  is  constrained  by  the  formula 
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b)  Ideally  shaped  asperities 


Fig.  6. 1 


Demonstration  of  modeling  technique  used  in  contact  area  formation 
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? N.  = AT  (6.1) 

J ' 


where  N = density  of  all  summit  penetrations  (equal  to  density  of  asperities) 

Nj  = density  of  summit  penetrations  of  5j 

If  the  corresponding  asperities  are  not  in  contact,  the  value  of  6 is  taken  to  be  negative. 

The  other  constraint  imposed  on  the  surfaces  in  contact  is  the  force  balance  that  must 
occur  at  the  interface.  According  to  the  selected  sign  for  6,  an  asperity  becomes  a load 
carrier  only  when  its  relative  penetration,  6,  is  positive.  In  general,  the  value  of  the 
contact  spot  radius  a and  the  load  F can  be  related  to  6 by  the  power  laws 

a.  = C,  , 6>0  (6.2) 

j i j i 

and 


F 

i 


0 


, 5.  < 0 

i 

,6.^0 

’ j 


(6.3) 


The  coefficients  Q are  dictated  by  the  assumed  asperity  shape  and  the  deformation 
mechanism.  Since  the  force  acting  at  the  interface  is  to  be  carried  by  the  discrete  contact 
spots,  one  can  write  the  force  relation  as 


E 

j 


N.F.  = 
i j 


P 


It  follows  that 


E 

j 


N.a  6 4 

; 3 j 


S;.k0 


= P 


(6.4) 
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To  perform  the  mathematical  procedures  analogous  to  those  in  Section  4.2,  the  two 
constraints,  Eqs.  6.1  and  4,  are  differentiated  as 

? dN.  = 0 (6.5) 

J 1 


and 


? dN.6c 4 

J 1 ' 


5y*0 


= 0 


(6.6) 


Suppose  that  there  exists  a statistical  weight  gj  which  represents  the  number  of  ways 
having  summit  penetration  6r  Application  of  Lagrange  multipliers  a and  p to  Eqs.  6.5 
and  6,  respectively,  leads  to  a distribution  law  as  follows 


Sje  a 

g.eae9&1 


, 6.  < 0 
’ j 

,6.^0 

’ j 


(6.7) 


The  variation  of  gj  is  contributed  by  some  interfacial  contact  features,  such  as  the 
surface  roughness  and  the  assumed  asperity  shape  and  deformation  manner.  For  that 
matter,  different  models  are  developed  to  model  an  interface  as  discussed  in  Section  3.2. 
Different  statistical  weights  are  proposed  as  formulated  in  the  next  four  sections. 

One  returns  to  the  Boltzmann  analysis.  The  two  multipliers  can  be  related  to  each 
other  by  substituting  Eq.  6.7  into  Constraint  1 (Eq.  6.1)  that  yields 


E 

j 


e ae 


0 


= N 
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so  that  e “ can  be  factored  out,  yielding 


N 


y i y.  p*,  i 

P + P p I \ 

j °j  I 6<  0 j Oj  l6. 


(6.8: 


Then,  substituting  e"“  into  Eqs.  6.7  and  4 gives  respectively 


N.  = 

j 


gjN 


ps,. 


S I + S 6 y'J>  I 

j °j  I 8y  < 0 y I 6^  2 0 


\7  -P  *,  ' 
gjNe  1 


y . y -p6c‘, 

Q + ' QP  J 

j bj  • 6 < 0 j &j  ' 6j*0 


,6.  < 0 
’ i 


,6.^0 

’ i 


(6.9) 


and 


^g.iVe155'  C 6.1'4|  n 

J °J  3 ; ISy20 

y:  i y,  p 6 c“ , 

j < o j °j  ' bj  * o 


(6.10) 


The  Boltzmann  statistical  model  in  conjunction  with  the  mechanical  analysis  has  been 
used  to  derive  the  distribution  of  the  summit  penetration  depths,  which  is  virtually  the 
distribution  of  the  contact  spots.  The  thermal  analysis  is  then  incorporated  to  obtain  the 
total  contact  conductance.  Including  the  effect  of  the  constriction  of  heat  flow  channels, 
given  by  Eqs.  3.1,  3.5,  and  6.2,  the  conductance  can  be  written  as 


h = — 


he 


iJj  J 


N.  C.  6 

i i i 


I 6^20 


(6.11) 
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It  is  assumed  that  all  contact  spots  are  evenly-spaced  such  that  the  ratio  of  the  spot  radius 
to  the  flow  channel  radius,  a/b  (see  Fig.  3.4),  is  the  same  for  each  spot.  Since  the 
constriction  alleviation  factor  Y is  a function  of  a/b,  Y can  be  factored  out  from  the 
summation  sign  as  shown  in  this  equation.  The  evaluation  of  the  ratio  a/b  will  be 
discussed  later  when  the  details  of  each  Boltzmann  statistics-based  model  are  described. 
If  Nj  in  Eq.  6.11  is  replaced  by  Eq.  6.9,  the  predicted  conductance  becomes 


h 


2k 


g.Ne'**'  C,6.qL 

J °J  l J 1 Sy2:0 


* E.i  . . 

j °J 1 Sy  < 0 j 6j  1 5y 


kO 


(6.12) 


When  all  the  discrete  variables  are  numerous  and  closely-spaced,  Euler-Maclaurin 
summation  formula  can  be  applied  to  change  the  summations  in  Eqs.  6.9,  10,  and  12  into 
integrals.  Then, 


dN, 

Ndb 


_ < 
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f S,d6  + f g6 
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(6.i3: 
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(6.14) 
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and 


Nf  g^C^db 


(6 . 15) 


The  three  equations  derived  above  are  one-to-one  functions  with  respect  to  the  Lagrange 
multiplier  p . It  can  thus  serve  as  an  intermediate  variable  which  forms  a link  among  the 
values  of  dN5/Nd6,  P,  and  h.  In  summary,  the  core  equations  for  a Boltzmann  statistics- 
based  model  are  formulated  and  they  can  be  solved  by  numerical  integration. 

As  discussed  earlier,  the  statistical  weight  function  g5  and  the  coefficients  Q depend 
on  the  surface  model  and  the  deformation  mechanism.  Since  there  are  several  possible 
ways  to  model  the  surface  roughness  and  the  asperity  deformation  (given  in  Section  3.2), 
g5  and  Q are  not  unique.  Four  different  models  are  analyzed  in  the  sections  that  follows. 


One  of  the  ways  to  represent  surface  irregularities  is  by  conical  asperities  and  the 
summits  of  the  touching  asperities  are  deformed  plastically.  Major  components  of  this 
modeling  technique  are  introduced  below. 

a)  Surface  model 

A rough  surface  is  modeled  by  a flat  plane  covered  with  cones  of  same  base  angle 


6.2  Plastic  Cone  Model 


and  different  heights  as  performed  by  Hsieh  [20];  see  Fig.  3.5.  Additionally,  assuming 
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a Gaussian  distribution  of  summit  heights  Hsieh  can  derive  a list  of  surface  characteristics 
relevant  to  this  model.  Some  of  the  characteristics  are  summarized  as  follows: 


1)  Density  of  surface  asperities  (number  per  unit  area), 


N = 


m 


7.308  a 


(6.i6: 


2)  Maximum  summit  height  (peak-to-valley  roughness), 


x =8  a 

max 


(6.17) 


3)  Mean  summit  height, 

[i  = 4a  (6.18) 

It  is  noted  that  the  asperities  having  heights  beyond  xmax  above  are  within  0.02%  of  the 
total  population  of  the  asperities  described  by  the  derived  Boltzmann  distribution  given  in 
Eq.  5.29. 

b)  Deformation  mechanism 

In  this  model,  the  asperities  are  deformed  following  the  Yovanovich  model  as 
illustrated  in  Fig.  6.2.  When  two  asperities  are  brought  into  contact,  their  summits 
(overlapping  regions)  penetrate  into  each  other  and  displace  plastically.  Similar  to  the  case 
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Fig.  6.2 


Plastic  deformation  at  asperity  summits 
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where  an  indenter  is  forced  into  a material  in  a hardness  test,  the  load  is  related  to  the  real 
contact  area  by 

F = HA  = Hna 2 (6.19) 

r 


From  the  geometry  given  in  Fig.  6.2,  the  contact  spot  radius,  a,  can  be  written  in 
terms  of  the  summit  penetration,  6 , as 

a = (6.20) 

2m 


so  that  Eq.  6.19  becomes 


F = H n 


2m 


7t  H 5.2 


4 m: 


(6.2i: 


Comparing  Eqs.  6.20  and  21  to  Eqs.  6.2  and  3,  respectively,  one  can  determine  the 


coefficients  as 


C, 

c2 

c3 


1 

2m 

1 

7 iH 
4 m2 


(6.22) 


C4  = 2 

c)  Statistical  weight 

The  product  of  the  statistical  weight  function  gs  and  d&  is  directly  proportional  to  the 
probability  of  having  a summit  penetration  depth  between  6 and  6 +d6.  Thus,  g5d6  can 
be  written  as 


ghd 6 = C x P( 6 to  6 +db) 


(6.23) 
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This  probability  involves  a combination  of  three  events.  They  are  the  height  of  an  asperity 
Xj,  the  height  of  the  opposing  asperity  x2,  and  the  base  plane  separation  they  are  related 
to  each  other  by 

6 = xx  + x2  - l (6.24) 

Thus, 

P(b  to  b+db)  = P{xx  fl  x2  fl  xx,x2  \ 'ito'i+d'i) 

Since  the  load  carrying  asperities  often  meet  at  their  summits  as  shown  in  Figs.  6.1 
and  2,  it  is  assumed  that  the  height  of  an  asperity  is  always  larger  than  its  deformed  region 
so  that  the  probability  can  be  more  specifically  expressed  by 

P(b  to  b+db)  = flx2^—  fl  jCj ,jc2 | ^to^+dE,) 

2 2 

For  any  given  xt  and  x2,  the  geometry  in  Fig.  6.2  shows  that  d6/d£  is  a constant  equal  to 
-1,  which  is  independent  of  6.  The  three  events  are  thus  mutually  exclusive  as 

P(b  to  b+db ) = P(x^~  n X2z-  fl  ltol+dl) 

which  is 

P(b  to  b+db)  = P(JC,*-)  x P(x2z-)  x P&tol+dl) 

2 2 


(6.25) 
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The  first  term  P(x,  * 6/2)  can  be  defined  as  the  portion  of  the  apparent  contact  area 
covered  by  asperities  of  height  greater  than  6/2.  By  using  the  probability  density  function 
given  in  Eq.  5.29,  the  following  can  be  written 


The  upper  bound  of  the  integral  is  the  maximum  summit  height,  8o,  which  is  given  by  Eq. 
6.17.  When  the  two  contacting  surfaces  have  identical  roughness,  P(x2  ;>  6/2)  equals 
P(x,  * 6/2). 

Owing  to  the  maximum  summit  height,  8o,  a contact  is  formed  only  if  the  base  plane 
separation  l is  smaller  than  16o,  which  serves  as  the  operating  domain  for  this  model. 
Thus,  the  last  term  is 


(6.26) 


pa  to  z+di)  - 


Moreover,  as  mentioned  earlier  that  d6/d?  is  -1,  the  above  equation  becomes 


pa  to  z.dt)  = - 


(6.27) 


Combining  Eqs.  6.23,  25,  26,  and  27  gives 
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As  given  in  Eqs.  6.14  and  15,  only  the  proportion  of  g5  is  of  interest.  Thus,  the  constant 
preceding  the  parenthesis  can  be  dropped,  which  will  not  affect  the  results  of  the  total 
conductance  versus  the  total  load.  It  is  thus  simplified  as 


This  expression  is  only  applicable  for  the  case  in  which  6 is  greater  than  zero, 
meaning  that  a contact  is  made.  When  6 is  less  than  zero,  gs  is  still  dependent  upon  the 
heights  of  the  asperities  involved  and  the  base  plane  separation.  Figure  6.3  illustrates  this 
situation.  Following  the  geometrical  relationship  among  xl5  x2,  l,  and  6 and  the  above 
procedures  deriving  g6,  one  can  show  that  1)  when  -8o<;6  <0, 


/ 


\ 2 


8o 


(6.28) 


\ 8/2 


/ 


(6.29) 


o 


o 


and  2)  when  -16o^6  <-8o, 


(6.30) 


o 


o 


The  statistical  weight  g6  is  thus  derived  for  the  entire  domain,  -16a<;6<:16o. 


d)  Constriction  alleviation  factor 


The  Gibson’s  expression,  listed  in  Table  3.1,  is  selected  to  determine  the  constriction 
alleviation  factor,  Y,  in  Eq.  6.15.  The  independent  variable  of  this  expression  is  the  ratio 
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Asperity  1 


Fig.  6.3 


Illustration  of  penetration  depth  for  non-contacting  asperities 
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of  contact  spot  radius  to  flow  channel  radius,  a/b,  and  it  can  be  determined  easily  for  this 
model.  In  this  effort,  one  can  denote  each  contact  spot  by  an  index  integer  i.  The  force 
balance  requires  that 


P_ 

H 


It  is  assumed  that  the  contact  spots  are  evenly  distributed  so  that  the  ratios  a/bj  are 
identical  for  all  spots.  Then, 


a_  P 
b \ H 


(6.31) 


In  this  equation,  the  loading  pressure  P has  been  given  in  Eq.  6.14. 
e)  Summary 

At  this  point,  all  parameters  involved  are  derived.  For  clarity  sake,  the  final  forms 
of  the  important  equations  are  summarized.  Substituting  coefficients  C,  in  Eq.  6.22  into 
Eqs.  6.14  and  15  gives 


and 


P = 


h = 


7 iH 
Am2 


k 

Wl|f 


16  o 

N J g&e  52dd 
o 

0 16o 

f * f g6e-*6ld 6 

-16o  0 

16o 

N [g&e  ^6d6 

0 

0 16o 

f Si,d 6 + f g6e'p&2d8 

-16a  0 


(6.32) 


(6.33) 
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where 


\ 

p 

H ) 

Combining  Eqs.  6.28,  29,  and  30  yields 


St  - ' 


(go  \ 2 
o o 


ff(x)x2dx  ,62:0 

V i/2 

160.6-*, 

ffCx^x2  J"  f(x2)x22dx2dxl  , -8 o £ 5 < 0 


o 

16o.6 


0 

16o.fi-*, 


J f(x j)x2  J f(x2)x22dx2dx1  , -16o^6<-8o 


(6.34) 


where 


1 / x - u Y 


f(x) 


2 \ a 


sjlft  o 


These  equations  conclude  the  development  of  the  Boltzmann  statistics-based  plastic  cone 
model. 


6.3  Plastic  Cone  and  Elastic  Base  Model 

Similar  to  the  plastic  cone  model  above,  the  roughness  of  a surface  can  be  modeled 
as  cones  with  their  summits  deformed  plastically.  In  addition,  based  on  Hsieh,  it  is 
assumed  that  the  base  of  a compressed  asperity  is  deformed  elastically  by  a deflection  r) 
as  illustrated  in  Fig.  6.4  [20].  In  this  figure,  the  dashed  lines  represent  the  original  shapes 
of  the  asperities  prior  to  deformation. 


67 


Fig.  6.4 


Illustration  of  plastic  summit  and  elastic  base  deformation 
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The  majority  of  the  equations  formulated  in  the  previous  section  remain  valid,  except 
d6/d£  which  is  no  longer  equal  to  -1.  The  geometry  given  in  Fig  6.4  shows  that 

6 = xl  + x2  - l - q,  - t]2  (6.35) 


This  difference  affects  the  expression  of  g6  for  6^0.  The  relationship  between  d£  and  d6 
must  be  re-established. 

The  elastic  base  deflection,  n,  can  be  approximated  by  the  result  of  a uniform 
distribution  of  pressure  acting  on  the  base  plane  over  the  asperity  base  area.  The  average 
deflection,  n , is  then  derived  and  given  by  [26] 


Using  Eq.  6.21, 


0.54(1  -\>2)mF 
Ex 


(6.36) 


T1  = 


0.54  rc  (1  - v2)Hb7 
AEmx 


(6.37) 


Thus,  l can  be  written  in  terms  of  6 by  substituting  Eq.  6.37  into  Eq.  6.35,  yielding 


\ = - 6 + xx  + x2  - 


0.547i//(l-v2)62  1 1 N 


AEm 


V *2  , 


(6.38) 


Differentiating  this  equation  with  respect  to  6 gives 


dl 

db 


0.27  nH(l-v2)b 
Em 


(6.39) 


69 


The  negative  sign  can  be  omitted  because  only  the  proportion  of  gs  is  of  interest.  Thus, 


dl  + 0.27  7tH(1-v2)S 


db 


Em 


V *1  + *2  / 


(6.40) 


Unlike  the  previous  model,  d6/d£  is  not  a constant  for  6^0.  The  statistical  weight 
function  cannot  be  simplified  and 


8o  8o 


S,  = ff(xi)  x,2  f f(x2 ) x2  dx2dxY 


5/2  5/2 


(6.41) 


For  the  entire  range  of  6 , g6  can  be  expressed  by 


S 5 


8 a 8a 


ffixjx2  ff(x2)x22-j?-dx2dx1  , 5z0 


6/2  6/2 


8o 


l6o  + 6-x 


jf(x j)Xj2  J f(x2)x22dx2dxl  , -8o<6<0 
0 0 

I60.6  16o*6-j:1 

J f(xx)x2  J f{x2)x2  dx2dxx  , -16o^6<-8o 


(6.42) 


In  summary,  the  conclusive  equations  for  this  model  include  Eqs.  6.32,  33,  and  42 


combined  with  Eq.  6.40  for  d£/d6 . 
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6.4  Elastic  Spherical  Asperity  Model 

Conical  asperities  have  been  employed  to  represent  the  irregularities  of  a rough 
surface.  Indeed,  several  investigations  demonstrate  that  modeling  the  roughness  by 
spherical  asperities  also  leads  to  reasonable  analytical  results.  One  example  is  the 
modified  Greenwood  and  Williamson  model,  developed  by  McWaid  and  Marschall  [24] 
and  briefly  discussed  in  Section  3.2.4.  In  the  present  Boltzmann  statistics-based  model, 


a great  deal  of  attributes  used  by  McWaid  and  Marschall  are  employed, 
discussed  below. 

Details  are 

a)  Surface  model 

The  techniques  of  random  process  theory  were  first  used  in  the  studies  of  ocean 
surfaces.  Nayak  then  extended  this  approach  to  analyze  the  statistics  of  a rough  surface 
by  using  the  power  spectral  density  of  a surface  profile  [36] . Comprehensive  derivations 
of  surface  geometric  parameters  can  be  found  in  the  publications  by  Nayak  [36]  and  Bush 
et  al.  [37],  Some  of  these  parameters  suitable  for  a predicting  model  for  contact 
conductance  are  the  density  of  asperities  (number  per  unit  area), 


mo. 

D = 4 

6 v/3  tc  mo2 

(6.43) 

the  mean  radius  of  summits,  R = 

8 sjmo. 

(6.44) 

the  mean  radius  of  summits, 


(6.44) 


71 


and  the  standard  deviation  of  summit  heights, 


0.8968 


P t 


\ 


mo, 


0 


(4.45) 


The  fundamental  variables  mo0,  mo2,  mo4,  and  p,  mentioned  before  in  Section  3.2.4, 
are  the  zeroth,  second,  fourth  moments  of  the  power  spectral  density,  and  the  bandwidth 
parameter,  respectively.  The  physical  meanings  of  mo0,  mo2,  and  mo4  are  the  variances 
of  the  surface  height,  the  slope,  and  the  curvature  of  a profilometric  sketch,  respectively. 

As  described,  mo0  is  the  square  of  the  standard  deviation  of  profile  height,  which  is 
commonly  obtained  by  a surface  profilometer.  Devices  measuring  the  root  mean  square 
of  surface  slope,  which  is  mo2'A,  although  not  as  popular,  are  also  commercially  available. 
The  major  issue  is  to  obtain  mo4  for  a surface.  The  practical  difficulties  of  this 
measurement  involve  the  choices  of  sample  spacing  and  signal  processing  [38], 

Regardless  of  the  stated  technical  difficulties,  mo2  and  mo4  can  still  be  obtained  by 
primitive  techniques.  Nayak  shows  that  not  only  the  density  of  asperities,  D,  but  also  the 
density  of  zero-crossings  Dzero  and  the  density  of  extrema  Dext  (peaks  and  valleys)  of  a 
profile  can  be  expressed  in  terms  of  mo2  and  mo4  [36],  Inversely,  mo2  and  mo4  are  written 
in  terms  of  Dzero  and  Dext  as 


mo.  = n 2monD  2 

2 0 zero 


(4.46) 


and 


(4.47) 
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Then,  by  counting  the  number  of  zero-crossings,  peaks,  and  valleys  of  a surface  profile, 
one  can  determine  mo2  and  mo4. 


b)  Deformation  mechanism 

When  two  spherical  asperities  are  pressed  together  and  meet  at  their  summits, 
Hertzian  theory  provides  a thorough  analysis  of  the  asperities  under  elastic  deformation 
[26],  Figure  6.5  depicts  the  deformation  scheme  of  this  kind.  The  original  shape  of  each 
asperity  is  drawn  by  the  dashed  line  and  the  elastically  deformed  asperities  are  drawn  by 
the  solid  lines.  For  two  contacting  spheres  having  the  same  radius,  Hertzian  theory 
predicts  that 

, 2 a2 


4Ea3 

3R(l  - v2) 


Recasting  the  6 equation  as  to 


N 2 


^ £0.5 


(6.48! 


and  using  it  in  the  F equation  above  gives 


J2RE 


. 1.5 


3(1  - v2) 


(6.49) 
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Fig.  6.5 


Material  deformation  between  spherical  asperities 
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From  these  two  equations,  the  coefficients  in  Eqs.  6.2  and  3 can  be  determined  as 


C2  = 0.5 


(6 . 50) 

c 4ire 

3 " 3(1  - v2) 

C4  = 1.5 


c)  Statistical  weight 

Similar  to  the  plastic  cone  model,  the  summit  heights,  Xj  and  x2,  the  base  plane 
separation,  £,  and  the  summit  separation,  6,  are  related  by  Eq.  6.24 

6 = x{  + x2  - ^ 

Therefore,  for  given  xl  and  x2,  d6/d£  equals  -1.  Following  the  reasoning  given  in  the 
plastic  cone  model,  Section  6.2,  the  statistical  weight  function  expressed  by  Eq.  6.34  can 
be  again  derived. 

d)  Constriction  alleviation  factor 

Since  it  is  assumed  that  for  all  contacts  the  spot  radius  to  flow  channel  radius  ratios 
are  equal  to  each  other,  this  leads  to 
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The  right  hand  side  of  this  equation  physically  means  the  square  root  of  the  ratio  of  true 
contact  area  to  apparent  contact  area.  Therefore,  a/b  can  be  alternatively  written  as 


Hence, 


(6.51) 


e)  Summary 

Substituting  the  coefficients  from  Eq.  6.50  into  Eqs.  6.14,  6.15,  and  6.51  results  in 
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(6.53) 
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and 
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T 
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N I g&e  ps'58d8 
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0 16  o 

d 6 + J gbe~?6'sd  8 


(6.54) 


f 8,ao  ♦ i g5e 


-16  o 


where  g6  is  obtained  by  Eq.  6.34  and  ¥(a/b)  is  obtained  by  the  Gibson's  expression  (see 
Table  3.1). 


6.5  Elastic  and  Plastic  Spherical  Asperity  Model 

In  addition  to  elastic  deformation,  continuous  loading  of  a spherical  asperity  may 
cause  a transformation  from  elastic  to  plastic  deformation  at  the  summit.  From  the  work 
of  Tabor  based  on  his  ball  indentation  hardness  test,  the  plastic  deformation  begins  when 
the  maximum  Hertzian  pressure  p0  reaches  60%  of  the  indentation  hardness  [27], 

p0  = 0.6  H (6.55) 

The  pressure,  p0,  occurs  at  the  center  of  the  contact  spot  and  its  value  relates  to  the  loading 
force,  F,  and  the  spot  radius,  a,  by  [26] 

3 F 

Po  - 


2na2 
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This  equation  can  be  expressed  in  terms  of  6 by  substituting  for  a and  F from  Eqs.  6.48 
and  49,  respectively.  Thus, 


J26E 

P0  = “ — 7~l=~ 

71(1  - \ )\[R 

Using  Eq.  6.55,  the  onset  of  plastic  deformation  becomes 


(6.56) 


s[2hE 


7X  (1  - \2)\[R 


= 0.6  H 


which  can  be  recast  as 


6 . 1.7765  »2(1  - ^ * 

E1 

The  plastically  deformed  region,  surrounded  by  elastic  material,  grows  as  the  summit 
penetration,  6,  increases  from  the  amount  given  above.  For  convenience,  a higher 
criterion  of  6 for  detectable  plastic  flow  is  selected  as  a critical  transition  from  a wholly 
elastic  mode  to  a wholly  plastic  mode  of  material  deformation.  Thus, 


. 2H2(1  - v2)2R 

8 = — (6.57) 

c E2 

The  summits  of  a pair  of  contacting  asperities  thus  undergo  elastic  or  plastic  deformation 
depending  on  their  penetration  depth,  6.  When  6 is  less  then  the  6C,  the  load  supported 
by  the  elastic  material  is  given  by  Eq.  6.49.  Once  6 exceeds  6C,  plastic  deformation 
becomes  dominant  so  that 


t i HR 


F = Hna2 
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(6.58) 
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In  summary,  the  loading  force  is 
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(6.59) 


Comparing  Eq.  6.59  to  Eq.  6.3,  one  can  determine  the  values  of  the  coefficients  C3  and 
C4  and  derive  them  in  two  groups  of  6 as  follows: 

1)  For  0^6  < 6C,  C3  and  C4  can  be  found  in  Eq.  6.50 


sflRE 
3(1  - v2) 
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2)  For  6*6C, 
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(6.60) 


Finally,  the  loading  pressure,  P,  and  the  total  thermal  contact  conductance,  h,  can  be 
written  by  substituting  all  the  above  coefficients  into  Eqs.  6.14  and  15  as 
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(6.61) 
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and 
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- 16o  0 «c 


(6.62) 


The  value  of  a/b,  needed  for  the  Gibson's  constriction  alleviation  factor  equation,  is 
obtained  from  Eq.  6.51  as 


a 


b 


"c  loo 

N Jg6e  ^6'5 6db  + IV  J g&e  p56d6 
■kR  o «c 

f g6d 8 + [g6e^6'5d6  + f g&e  Pt>d5 

-16o  0 6C 


(6.63) 


The  statistical  weight  g6  remains  the  same  as  in  Eq.  6.34. 


CHAPTER  7 

NUMERICAL  RESULTS  AND  DISCUSSION 

The  equations  necessary  to  predict  the  thermal  contact  conductance  have  been 
formulated  in  the  previous  chapter.  Finding  the  analytical  solutions  to  these  equations  is 
not  successful  because  of  the  difficulty  of  evaluating  the  multiple  integrations  involved. 
For  this  reason,  numerical  methods  are  used.  In  the  following,  the  computation  technique 
and  the  numerical  results  are  presented  and  discussed. 

7.1  Numerical  Method 

Essentially  all  definite  integrals  are  evaluated  numerically  by  Simpson's  rule.  The 
only  exception  is  the  statistical  weight  g&  given  by  Eq.  6.28  which  can  be  simplified 
analytically  to  an  expression  in  terms  of  the  error  function.  The  relevant  equation  for  this 
derivation  is  Eq.  A. 3,  provided  in  Appendix  A. 

When  Simpson's  rule  is  applied,  an  initial  and  arbitrary  number  of  equally-spaced 
increments  of  the  dummy  variable  of  integration  are  selected  to  perform  a numerical 
integration.  Then,  the  computation  is  repeated  by  using  twice  the  number  of  increments. 
In  other  words,  the  size  of  each  increment  is  reduced  by  one  half.  This  procedure  is 
repeated  until  the  termination  criterion  is  reached.  The  termination  criterion  applied  here 
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is  that  the  errors  between  consecutive  computations  do  not  exceed  one  percent.  It  is  found 
that  a sufficient  number  of  increments  required  to  produce  a satisfactory  solution  varies 
from  100  to  400. 


7.2  Numerical  Data 


Experimental  data  from  the  open  literature  have  been  used  for  comparison  with  the 
numerical  results  predicted  by  the  Boltzmann  statistics-based  models.  Only  those 
experiments  conducted  in  a vacuum  environment  with  the  necessary  topographic 
parameters  measured  and  reported  are  included  for  comparison.  The  surface 
characteristics  and  the  sources  of  reference  for  all  the  selected  data  are  summarized  in 
Table  7.1. 

Since  the  surface  slope  is  not  commonly  reported  in  the  literature,  only  eleven  sets 
of  qualifying  data  are  found  useful  and  they  cover  a wide  range  of  surface  roughness  and 
material  properties.  The  specimens  include  Stainless  Steels  303,  304,  and  416,  Nickel 
200,  Zirconium-2. 5 wt%  Nb,  and  Zircaloy-4  and  the  surface  roughness  varies  from  0.338 
pun  (smooth)  to  9.86  pirn  (rough).  The  results  are  plotted  in  Figs.  7.1  through  11. 

Useful  experimental  data  to  assess  the  accuracy  of  the  spherical  models  is  even 
scarcer.  Measurements  of  both  the  thermal  contact  conductance  and  the  moments  of 
power  spectral  density  are  found  only  in  the  study  by  McWaid  and  Marschall  [24],  The 
specimens  tested  were  made  of  Aluminum  alloy  6061  and  Stainless  Steel  304  with  two 
types  of  finishes:  bead-blasted  and  ground.  The  four  sets  of  data  provided  by  this 


Table  7. 1 Experiments  for  thermal  contact  conductance  of  nominally-flat  surfaces  in  a vacuum  environment 
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Fig.  7.1 


Thermal  contact  conductance  versus  contact  pressure  for  Stainless  Steel  303 
surfaces  in  a vacuum  (Data  1) 


Contact  Conductance,  h (Btu/  hr  F ft  ) 
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Fig.  7.2  Thermal  contact  conductance  versus  contact  pressure  for  Stainless  Steel  303 
surfaces  in  a vacuum  (Data  2) 


Contact  Conductance,  h (Btu/  hr  F ft  ) 
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Fig.  7.3  Thermal  contact  conductance  versus  contact  pressure  for  Stainless  Steel  416 
surfaces  in  a vacuum  (Data  3) 


Contact  Conductance,  h (W/m  K) 
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Fig.  7.4  Thermal  contact  conductance  versus  contact  pressure  for  Nickel  200  surfaces 
in  a vacuum  (Data  4) 


Contact  Conductance,  h (W/m  K) 
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Fig.  7.5  Thermal  contact  conductance  versus  contact  pressure  for  Nickel  200  surfaces 
in  a vacuum  (Data  5) 


Contact  Conductance,  h (W/m  K) 
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Fig.  7.6  Thermal  contact  conductance  versus  contact  pressure  for  Stainless  Steel  304 
surfaces  in  a vacuum  (Data  6) 


Contact  Conductance,  h (W/m  K) 


89 


Fig.  7.7  Thermal  contact  conductance  versus  contact  pressure  for  Stainless  Steel  304 
surfaces  in  a vacuum  (Data  7) 


Contact  Conductance,  h (W/m  K) 
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Fig.  7.8  Thermal  contact  conductance  versus  contact  pressure  for  Zirconium- 
2.5  %wt  Nb  surfaces  in  a vacuum  (Data  8) 


Contact  Conductance,  h (W/m  K) 
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Fig.  7.9  Thermal  contact  conductance  versus  contact  pressure  for  Zirconium- 
2.5  %wt  Nb  surfaces  in  a vacuum  (Data  9) 


Contact  Conductance,  h (W/m  K) 
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Fig.  7. 10  Thermal  contact  conductance  versus  contact  pressure  for  Zircaloy-4  surfaces 
in  a vacuum  (Data  10) 


Contact  Conductance,  h (W/m  K) 
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Fig.  7.11  Thermal  contact  conductance  versus  contact  pressure  for  Zircaloy-4  surfaces 
in  a vacuum  (Data  11) 
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reference  are  used  to  validate  the  spherical  models.  The  test  results  of  the  bead-blasted 
surfaces  are  presented  in  Figs.  7.12  and  13.  Tests  of  the  contact  between  the  ground 
surfaces  oriented  with  lays  parallel  and  perpendicular  to  each  other  are  plotted  in  Figs. 
7.14  and  15,  respectively.  It  should  be  noted  that  these  figures  (Figs.  7.1  through  15)  also 
include  the  predictions  made  by  others  as  reported  in  the  literature  (Section  3.2). 

A dimensionless  plot,  Fig.  7.16,  using  parameters  derived  in  the  previous  correlation 
studies  [19,29],  further  compares  the  prediction  by  the  plastic  cone  model  and  the  first 
eleven  sets  of  the  experimental  data  listed  in  Table  7.1.  In  this  comparison,  the  data  are 
plotted  in  dimensionless  contact  conductance, 

, . ha 
h = 

mk 


versus  the  dimensionless  pressure, 


P-  = 


_P 

H 


Elastic  deformation  is  accounted  for  in  the  other  three  Boltzmann  statistics-based  models 
(plastic  cone  and  elastic  base,  elastic  spherical  asperity,  elastic  and  plastic  spherical 
asperity).  This  leads  to  the  presence  of  the  modulus  of  elasticity  and  the  Poisson  ratio  in 
the  analysis.  A dimensionless  plot  becomes  inappropriate  in  such  models. 


Contact  Conductance,  h (W/m  K) 
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Fig.  7. 12  Thermal  contact  conductance  versus  contact  pressure  for  Aluminum  6061-T6 
surfaces  in  a vacuum  (Data  12) 


Contact  Conductance,  h (W/m  K) 
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Fig.  7.13  Thermal  contact  conductance  versus  contact  pressure  for  Stainless  Steel  304 
surfaces  in  a vacuum  (Data  13) 


Contact  Conductance,  h (W/m  K) 
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Fig.  7. 14  Thermal  contact  conductance  versus  contact  pressure  for  Stainless  Steel  304 
surfaces  in  a vacuum  (Data  14) 


Contact  Conductance,  h (W/m  K) 
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Fig.  7.15  Thermal  contact  conductance  versus  contact  pressure  for  Stainless  Steel  304 
surfaces  in  a vacuum  (Data  15) 
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P* 


Fig.  7.16 


Dimensionless  thermal  contact  conductance  versus  dimensionless  contact 
pressure 
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7.3  Discussion 

The  numerical  results  by  the  conical  models  are  compared  to  the  contact  conductance 
data  measured  by  Mikic  and  Rohsenow  [39]  and  Henry  [40]  in  Figs.  7.1,  2,  and  3.  The 
difference  between  the  predicted  and  the  empirical  values  is  within  30%.  The  predictions 
by  the  Hsieh  model  and  the  Yovanovich  model  are  also  presented  in  these  three  figures. 
In  Fig.  7.1,  the  Yovanovich  model  overpredicts  the  contact  conductance  while  the  Hsieh 
model  and  the  conical  models  are  in  good  agreement  with  the  experimental  data. 
However,  in  Figs.  7.2  and  3,  the  conical  models  are  found  to  be  the  best  among  them. 

The  experiments  conducted  by  Hegazy  [31]  are  presented  in  Figs.  7.4  through  11. 
Here,  the  trend  observed  is  consistent  in  that  the  predictions  by  the  conical  models  agree 
well  with  the  experiment  at  low  loading  pressures.  As  the  pressure  increases,  the  models 
underpredict  the  contact  conductance.  A comparison  with  the  Yovanovich  model  indicates 
that  the  conical  models  work  better  at  low  loading  pressures  while  the  Yovanovich  model 
works  better  at  high  loading  pressures. 

It  is  found  that  the  predicted  values  of  the  thermal  contact  conductance  by  the  plastic 
cone  and  elastic  base  model  are  consistently  higher  than  those  by  the  plastic  cone  model. 
As  the  two  predictions  are  very  close,  some  numerical  data  are  plotted  again  with  an 
enlarged  scale  in  Fig.  7.17  to  illustrate  the  claim.  This  phenomenon  is  explained  by  the 
fact  that  the  elastic  base  deflection  of  a load  carrying  asperity  is  considered  in  the  plastic 
cone  and  elastic  base  model.  Thus,  for  a given  contact  pressure,  even  though  the  total 
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Contact  Pressure,  P (Pa) 


a)  Data  5 (E/H  = 95.4) 


Fig.  7.17  Thermal  contact  conductance  versus  contact  pressure  for  a)  Data  5 and 
b)  Data  6 


Contact  Conductance,  h (W/m  K) 
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b)  Data  6 (E/H  = 50.4) 


Fig.  7. 17  — continued 


103 


actual  contact  area  is  the  same,  more  contact  spots  are  formed  as  illustrated  in  Fig.  7.18. 
Hence,  the  resultant  conductance  is  higher. 

One  expects  that  if  the  modulus  of  elasticity,  E,  is  large,  its  effect  on  the  prediction 
of  contact  conductance  would  be  small  as  a results  of  the  small  base  deflection.  This 
variation  can  be  found  in  the  numerical  results.  Among  all  the  data  presented  in  Fig.  7.1 
through  Fig.  7.7,  the  two  extreme  values  of  the  ratio  E/H  are  95.4  and  50.4  for  Data  5 
and  Data  6,  respectively.  The  plots  in  Figs.  7.17a  and  17b,  set  up  in  the  same  scale, 
indicate  that  the  smaller  the  E/H  is,  the  larger  the  difference  between  the  two  predictions 
is.  Thus,  the  numerical  results  agree  with  the  physical  expectation. 

In  Figs.  7.12  and  13,  the  predictions  by  the  spherical  models  are  compared  to  the 
empirical  data  obtained  from  the  tests  of  bead-blasted  aluminum  and  stainless  steel 
surfaces.  The  disagreement  is  found  to  be  large.  However,  the  modified  Greenwood  and 
Williamson  model,  also  using  spherical  asperities,  appears  to  be  satisfactory  in  predicting 
the  contact  conductance. 

The  comparisons  for  the  ground  surfaces  are  presented  in  Figs.  7.14  and  15.  As  the 
surfaces  are  strongly  anisotropic,  the  equivalent  values  of  the  spectral  moments  given  by 
McWaid  and  Marschall  [24]  are  used  in  the  prediction.  Again,  the  present  spherical 
models  overpredict  the  contact  conductance  in  all  cases  tested. 

The  elastic  asperity  model  tends  to  predict  a lower  conductance  than  the  elastic  and 
plastic  asperity  model  at  high  contact  pressure  as  shown  in  Figs.  7.13  to  15.  This  can  be 
ascribed  to  the  fact  that  at  high  pressure  many  touching  asperities  should  be  deformed 
plastically  while  the  elastic  spherical  asperity  model  still  analyzes  the  contact  by 
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a)  Plastic  cone  model 


b)  Plastic  cone  and  elastic  base  model 


Sane  total  contact  area 
but  nore  contact  spots 


Fig.  7.18  Illustration  of  the  difference  in  deformation  of  asperities  for  a given 
contact  pressure  between  a)  the  plastic  cone  model  and  b)  the  plastic  cone 
and  elastic  base  model 
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considering  the  asperities  in  elastic  deformation.  The  stress-strain  diagrams  given  in  Fig. 
7.19  provide  the  basis  that  less  contact  area  is  predicted  by  the  elastic  spherical  asperity 
model  at  high  contact  pressure.  As  a result,  the  estimates  of  the  thermal  contact 
conductance  are  low. 

The  dimensionless  plot  given  in  Fig.  7.16  effectively  summarizes  the  data  discussed 
in  this  chapter.  Here,  all  the  predictions  by  the  plastic  cone  model  are  plotted.  The 
correlation  between  the  predicted  IT  and  P*  is  expressed  by 

h‘  = 0.170  P*0  693  (7.1) 

The  good  agreement  (±30%)  between  the  prediction  by  the  plastic  cone  model  and 
the  measurements  by  Mikic  and  Rohsenow  [39]  and  Henry  [40]  spans  a wide  range  of  P* 
(4xl0‘4  < P*  < 3xl02).  For  Hegazy's  [31]  experiments,  good  agreement  (±30%)  is 
found  when  P*  is  less  than  6xl04.  When  P*  is  greater  than  6xl0  4,  the  model  tends  to 
underpredict  and  the  deviation  increases  with  P*.  The  worst  case  found  is  at  P* 
approximately  equal  to  4xl0'3,  where  the  experimental  h*  is  twice  as  high  as  the  predicted 
h*. 

On  the  basis  of  the  above  analysis,  the  predictions  of  thermal  contact  conductance  by 
the  conical  models  at  low  loads  (P*  < 6x1  O'4)  are  highly  reliable.  At  high  load,  since 
there  is  a large  deviation  among  existing  empirical  data,  it  is  difficult  to  assess  the 
accuracy  of  the  predicting  models.  Nonetheless,  the  predictions  by  the  conical  models  are 
within  the  uncertainty  in  the  experimental  data. 
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a)  Elastic  model 


b)  Elastic  and  plastic  model 


Fig.  7.19  Stress-strain  analyses  used  by  the  two  spherical  asperity  models: 
a)  elastic  model  and  b)  elastic  and  plastic  model 
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In  summary,  the  Boltzmann  statistical  model  is  found  to  perform  well  with  the 
assumptions  made  in  the  two  conical  models.  The  deviation  between  their  predictions  are 
negligibly  small  for  all  the  cases  tested.  In  addition,  the  consideration  of  the  elastic  base 
deformation  complicates  the  computation  and  that  costs  a longer  computer  time  to  obtain 
a numerical  solution.  Thus,  it  is  concluded  that  the  plastic  cone  model  is  the  most 
effective  among  all  the  Boltzmann  statistics-based  models  developed  in  this  investigation. 


CHAPTER  8 

MEASUREMENT  OF  ELECTRICAL  CONTACT  CONDUCTANCE 


Experiments  were  conducted  in  order  to  provide  additional  data  to  evaluate  the 
accuracy  of  the  Boltzmann  statistics-based  models  in  the  prediction  of  the  contact 
conductance.  However,  instead  of  measuring  the  thermal  contact  conductance,  which 
requires  a costly  apparatus,  the  electrical  contact  conductance  was  measured  as  the  two 
conductances  appear  to  be  analogous. 

Some  investigators  have  studied  the  electrical  contact  conductance  between  metals, 
such  as  silver,  aluminum,  copper,  low  carbon  steel,  iron,  among  others  [41,42],  It  has 
been  shown  that  the  surface  film  on  metal,  such  as  oxide  layer,  yields  unnoticeable  effect 
on  the  thermal  conductance  yet  reduces  the  electrical  conductance  considerably.  For  this 
reason,  stainless  steel  was  chosen  to  be  the  test  material  as  it  is  highly  resistant  to 
corrosion.  Moreover,  the  low  electrical  conductivity  of  the  stainless  steel  (approximately 
2x  106  Q'nr1)  is  advantageous  from  the  standpoint  that  the  minute  voltage  drop  across  the 
contact  interface  can  be  more  sensitively  measured  by  instruments.  The  details  of  the 
experimental  procedures  and  the  test  results  are  given  below. 
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8.1  Surface  Preparation  and  Measurement 

Two  sets  of  stainless  steel  specimens  (a  pair  in  each  set)  were  prepared  with  different 
surface  textures.  Each  specimen  was  cut  from  a 25  mm  diameter  Stainless  Steel  304  bar 
stock  and  then  machined  down  to  19  mm  long  by  milling  and  turning.  Its  surfaces  were 
lapped  with  a 600  grit  compound  to  reduce  waviness.  The  flatness  measurement  was  made 
by  using  a 20  mm  surface  profile  trace  drawn  by  a Dektak  profilometer.  It  was  found  that 
the  surface  was  originally  spherical  in  shape  with  the  center  being  higher.  The  waviness 
readings  from  the  peak  to  the  valley  were  greater  than  1 /zm.  In  order  to  reduce  this 
waviness,  the  lapping  technique  used  to  produce  gauge  blocks  was  applied.  In  this 
technique,  a set  of  three  surfaces  were  lapped  against  each  other  (two  at  a time)  alternately 
and  repeatedly  with  fine  lapping  compounds  (800  grit,  1000  grit,  and  1200  grit 
successively).  The  surfaces  were  measured  again  and  found  that  the  flatness  was  improved 
to  a waviness  height  of  less  than  0.7  /zm. 

Surfaces  were  roughened  by  sandblasting  with  beads  of  sizes  between  0.5  and  0.8 
mm.  The  specimen  was  placed  20  cm  away  from  the  nozzle  of  the  blasting  gun.  The 
desirable  roughness  was  accomplished  by  adjusting  the  blaster  pressure  (60  psi  for  rough 
surfaces  and  40  psi  for  smooth  surfaces).  Each  surface  was  sandblasted  until  a 
homogeneous  texture  was  found  under  a microscope.  The  total  time  required  for  this 
blasting  process  was  about  15  seconds. 

The  processed  surfaces  were  then  cleaned  ultrasonically  with  detergent  solution, 
acetone,  and  distilled  water  successively.  Again,  for  each  surface  two  20  mm  profiles 
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perpendicular  to  each  other  were  recorded  by  a Dektak  profilometer  (see  figures  in 
Appendix  B).  In  these  recorded  profiles,  the  waviness  to  roughness  ratio  is  found  to  be 
small  to  the  extent  that  the  surface  waviness  is  negligible.  The  root-mean-square 
roughness  of  the  surface  was  measured  by  a Brush  profilometer  and  their  values  are 
summarized  in  Table  8.1. 

8.2  Apparatus  Setup  and  Data  Collection 

The  material  properties  of  the  stainless  steel  bar  stock  had  been  measured  before  the 
electrical  contact  conductance  was  tested.  The  Vickers  hardness  number  was  found  to  be 
299  (425  kpsi)  by  a Buchler  Micromet  II  microhardness  tester.  The  electrical  conductivity 
of  the  steel  was  measured  by  a setup  illustrated  in  Fig.  8.1.  Here,  a 1.2-V,  40  amp  hour 
battery  was  used  to  produce  a 1 A current  during  the  test.  The  electrodes  were  attached 
to  the  stainless  steel  bar  by  means  of  a brass  ring  where  three  screws  are  provided  as 
shown  in  Fig.  8.2.  The  voltage  drop  across  the  electrodes  separated  by  a distance  L was 
of  the  order  of  10  \ aV,  which  is  too  small  for  conventional  voltmeters  to  measure. 
Therefore,  a Keithley  Model  182  digital  voltmeter  (0.01  /xV  resolution)  was  used.  The 
measured  resistance  Re  versus  the  separation,  L,  and  a linear  regression  line  are  plotted 
as  shown  in  Fig.  8.3.  The  slope  of  the  line  was  determined  to  be  1 .47  x 10'3  Qm4,  which 
yields  an  electrical  conductivity  of  1.34X106  Q ’m1.  The  linear  relationship  among  the 
data  points  also  indicates  the  reliability  of  this  techniques  in  measuring  low  resistance. 


Table  8.1  Summary  of  surface  roughness 


Root-mean-square  roughness 

low  (jum) 

high  (jim) 

average  (jum) 

Set  1 

Specimen  A 

1.372 

1.473 

1.423 

Specimen  B 

1.143 

1.473 

1.308 

Set  2 

Specimen  C 

0.8636 

1.016 

0.9398 

Specimen  D 

0.6604 

0.8128 

0.7366 
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Fig.  8.1 


Schematic  of  electrical  conductivity  measurement 
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Side  view  (Section  A-A) 


specimen 


Fig.  8.2 


Connection  of  an  electrode  to  a specimen 


(10  Q) 
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L (m) 


Fig.  8.3  Measured  resistance  of  stainless  steel  bar  stock  versus  length 
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A test  rig  was  designed  for  the  measurement  of  the  electrical  contact  conductance  as 
depicted  in  Fig.  8.4.  The  features  of  this  setup  are  analogous  to  those  of  the  schematic  for 
a thermal  test  shown  in  Fig.  3.9.  The  electrodes  were  attached  to  the  copper  bars  and  the 
stainless  steel  specimens  in  the  same  manner  as  before  (Fig.  8.2).  The  copper  was  chosen 
for  its  softness  and  high  electrical  conductivity.  It  allows  the  current  to  flow  uniformly 
through  the  specimens.  The  stack  was  placed  between  two  platforms  of  a hydraulic  press. 
At  the  top  end  of  the  stack,  the  test  pressure  was  monitored  by  means  of  a load  cell 
constructed  and  calibrated  by  the  Civil  Engineering  Department  of  the  College  of 
Engineering  of  the  University  of  Florida.  A ball  bearing  was  also  used  such  that  any 
misalignment  of  the  stack  would  not  disturb  the  loading  at  the  interface.  Finally, 
insulators  were  used  so  that  the  current  would  not  be  channeled  through  the  hydraulic 
press. 

During  the  test  the  rheostat  was  adjusted  to  maintain  a 1 A current  in  the  circuit. 
Every  time  the  load  was  changed,  a waiting  period  of  ten  minutes  was  scheduled  to  allow 
the  instrument  readings  to  stabilize.  Then,  the  voltage  drop  V across  the  contact  interface, 
the  current  I,  and  the  load  cell  signal  were  recorded. 

The  measured  resistance  Remeas,  equal  to  V/I,  includes  both  the  contact  resistance 
and  the  bulk  resistance  R,.  bulk.  Then,  the  contact  resistance  is  obtained  by  using  the 
formula 


where 


R 


e,cont 


R 


e,meas 


R 


e.bulk 


(8.1) 


e,meas 
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insulator 

hydraulic 
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Fig.  8.4 


Schematic  of  the  apparatus  for  testing  electrical  contact  conductance 
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and 


R 


e,bulk 


_L^ 

kA 


The  corresponding  contact  conductance  per  unit  area  is 


h 


e 


l 


R ,A 

e,cont 


(8.2) 


The  empirical  results  are  plotted  and  discussed  in  the  following  section. 

8.3  Empirical  Results  and  Discussion 

The  measured  values  of  the  electrical  contact  conductance  are  presented  in  Figs.  8.5 
and  6.  The  major  source  of  error  is  the  uncertainty  in  the  measurement  of  L,  which  is 
used  to  determine  Re  bulk.  As  he  is  related  to  Rcbulk  by  Eqs.  8.1  and  2,  the  magnitude  of  the 
error  increases  with  the  applied  load. 

The  experimental  results  are  compared  to  the  prediction  by  the  plastic  cone  model  as 
it  was  found  to  be  the  most  satisfactory  among  all  the  Boltzmann  statistics-based  models. 
It  is  found  that  the  measured  values  of  the  electrical  contact  conductance  are  considerably 
lower  than  the  predicted  values  by  Eq.  6.33.  A correction  factor  Xe  is  introduced  to 
improve  the  predicted  values  and  the  results  are  plotted  with  the  experimental  data.  The 
values  of  A,e  are  0.0302  and  0.0704  for  the  sets  1 and  2,  respectively. 


Electrical  Contact  Conductance, 
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Fig.  8.5  Electrical  contact  conductance  versus  contact  pressure  for  Stainless  Steel  304 
surfaces  (Set  1) 


Electrical  Contact  Conductance, 
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Fig.  8.6  Electrical  contact  conductance  versus  contact  pressure  for  Stainless  Steel  304 
surfaces  (Set  2) 
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One  can  observe  that  the  measured  electrical  contact  conductance  increases  as  the 
surface  roughness  increases.  This  phenomenon,  also  reported  in  references  [42,43],  is  in 
contrast  with  the  traditional  relationship  between  the  contact  conductance  and  the  surface 
roughness  because  of  the  presence  of  surface  films.  In  the  experiment  by  Tsukizoe  and 
Hisakado  [42],  both  electrical  and  thermal  contact  conductances  were  tested.  The  lowest 
and  highest  values  of  the  dimensionless  electrical-to-thermal  ratio,  equivalent  to 
presented  here,  for  a contact  between  copper  surfaces  was  found  to  be  0.0476  for  smooth 
surfaces  and  0.143  for  rough  surfaces,  respectively. 

The  physical  significance  of  A,e  is  related  to  the  effect  of  the  presence  of  surface  film, 
which  is  electrically  insulating.  In  general,  a metal  surface  is  covered  by  a composition 
of  tarnish  and  alien  films  [44],  A tarnish  film,  such  as  an  oxide  or  sulphide,  is  a film  in 
which  one  constituent  consists  of  the  atoms  of  the  substrate  material.  An  alien  film,  on 
the  other  hand,  is  formed  by  substances  in  the  surroundings  deposited  on  a metal  surface, 
for  example,  water  or  lubricant  film.  The  surface  films  permit  the  flow  of  phonons  but 
inhibit  the  flow  of  electrons.  The  reason  is  that  the  film  and  the  metal  have  similar 
thermal  conductivities  but  the  electrical  conductivity  of  the  film  is  much  lower  than  that 
of  the  metal.  For  instance,  the  electrical  conductivity  of  copper  oxide  is  5xl0'3  flr'm'1  and 
that  of  pure  copper  is  5 X 107  Q ’m1.  Thus,  when  two  surface  asperities  are  in  contact,  the 
electrons  can  be  conducted  through  a portion  of  the  contact  spot  area  only  where  the 
surface  films  are  ruptured.  This  results  in  a low  electrical  contact  conductance  compared 


to  the  thermal  one. 
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One  surface  film  rupture  mechanism  is  by  asperity  penetration  as  illustrated  in  Fig. 
8.7  [43].  This  figure  shows  that  under  a given  load  the  sharper  asperity  (found  in  a 
rougher  surface)  produces  a larger  pure  metal-metal  contact  area  (electrically  conductive) 
but  the  same  load-carrying  area  (thermally  conductive).  This  explains  why  the  contact 
conductance  decreases  with  the  surface  roughness  in  a thermal  test  but  increases  in  an 
electrical  test  as  the  values  of  Xe  are  found  in  Figs.  8.5  and  6. 

Another  mechanism  that  induces  pure  metal-metal  contact  is  by  film  cracking.  Osias 
and  Tripp  performed  a test  modeling  a film  covered  spherical  asperity  of  a metal  surface 
by  a plastic  sphere  coated  with  a brittle  shell  [45],  Subsequent  to  the  loading  of  this 
sphere,  the  modes  of  cracking  observed  are  radial,  circumferential,  and  craze  cracks. 
Thus,  in  reality  electrons  can  be  conducted  through  these  cracks  of  an  insulating  film. 

If  a surface  film  is  sufficiently  thin  (thickness  less  than  50  A),  electrons  can  pass 
through  the  film  by  tunneling  effect.  In  this  mechanism,  an  electron  having  a potential 
higher  than  the  work  function  of  the  film  will  jump  over  the  potential  barrier.  Thus,  a 
flow  of  electrons  is  induced. 

In  this  experimental  investigation,  it  is  found  that  the  surface  film  effect  on  the 
electrical  contact  conductance  between  two  stainless  steel  surfaces  is  significant.  The 
magnitude  of  the  correction  factor,  A,e,  is  comparable  to  those  found  in  the  literature  [42] . 


122 


a) 


b> 


Fig.  8.7 


Demonstration  of  larger  pure  metal -metal  contact  area  produced  by 
a)  rough  surfaces  in  contact  than  b)  smoother  surfaces  in  contact 


CHAPTER  9 

CONCLUSIONS  AND  RECOMMENDATIONS 


This  investigation  has  shown  that  the  prediction  of  thermal  contact  conductance  by 
statistical  mechanics  is  viable.  Considering  surface  asperities  comparable  to  molecules  in 
the  traditional  sense  of  microscopic  viewpoint,  the  Boltzmann  statistical  model  successfully 
determines  the  "most  probable  distribution"  of  the  contact  spots.  Thus,  a relationship 
between  the  contact  conductance  and  the  relative  pressure  is  derived. 

A thorough  study  was  conducted  by  applying  various  assumptions  of  the  surface 
model  and  the  deformation  mechanism  in  conjunction  with  the  computation  procedures  of 
the  Boltzmann  model.  It  is  found  that  modeling  a rough  surface  by  conical  asperities 
whose  summits  are  deformed  plastically  is  the  most  effective  approach.  The  disagreement 
between  the  predicted  thermal  conductance  and  the  published  data  is  less  than  30%  when 
the  dimensionless  pressure  (P/H)  is  less  than  6xl0'4.  Over  the  entire  test  domain,  the 
worst  empirical  data  is  twice  as  much  as  the  prediction. 

The  results  of  the  models  using  spherical  asperities  are  less  accurate.  However,  the 
prediction  by  the  elastic  and  plastic  spherical  asperity  model  gives  a trend  that  is 
comparable  to  the  experimental  data. 

The  differences  in  the  predictions  by  the  four  Boltzmann  statistics-based  models 
developed  agree  with  the  expectations.  For  instance,  due  to  the  base  deformation 
mechanism,  the  plastic  cone  and  elastic  base  model  predicts  a thermal  contact  conductance 
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that  is  higher  than  the  plastic  cone  model.  Moreover,  the  elastic  spherical  asperity  model 
predicts  a conductance  that  is  lower  than  the  elastic  and  plastic  asperity  model.  This  is 
because  the  former  overestimates  the  stress  on  a contact  under  high  load. 

A number  of  recommendations  are  suggested  to  explore  this  approach  further  in  the 
future.  First,  the  statistical  weight  function,  g6,  is  defined  on  the  basis  of  the  base  plane 
separation,  E,  varying  between  0 and  16o  (2  times  the  highest  asperity).  In  reality,  if  E, 
becomes  less  than  a critical  value  £c,  mechanical  failure  of  the  bulk  materials  occurs. 
Thus,  it  is  more  reasonable  that  \ varies  between  Ec  and  16o.  The  value  of  Ec  can  be 
determined  by  the  existing  models  that  relate  E,  and  the  loading  pressure.  Candidate 
models  include  the  Hsieh  model  [20]  and  the  modified  Greenwood  and  Williamson  model 
[24], 

The  Boltzmann  statistical  model  has  been  explored  in  this  investigation  and 
reasonably  accurate  predictions  are  obtained.  The  Bose-Einstein  and  Fermi-Dirac 
statistical  models,  assuming  that  the  particles  are  indistinguishable,  have  been  applied 
successfully  in  statistical  thermodynamics.  Thus,  these  two  models  may  also  play  an 
important  role  in  the  field  of  contact  conductance.  Investigations  of  these  two  statistical 
models  are  recommended. 

A correlation  among  the  numerical  results  of  the  contact  conductance  and  the 
parameters  involved  is  also  encouraged.  A successful  correlation  should  simplify  the 
complicated  numerical  integrations. 

The  experimental  results  of  this  research  show  that  the  presence  of  electrically 
insulating  surface  films  significantly  reduce  the  pure  metal-metal  contact  area  between  the 
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stainless  steel  surfaces.  The  film  effect  results  in  a low  measured  electrical  contact 
conductance.  Film  penetration  and  film  cracking  are  the  fundamental  mechanisms  that 
produce  pure  metal-metal  contacts,  where  the  electrons  are  easily  transferred.  The 
electrons  can  also  be  transferred  by  the  tunneling  effect  if  the  film  is  sufficiently  thin.  At 
the  present  time,  the  models  that  predict  the  film  effect  partially  depend  on  experimental 
measurements. 

Still,  testing  contact  conductance  electrically  plays  a vital  role  in  the  field  of  thermal 
contact  conductance.  The  main  reason  is  the  cost  of  a test  facility  being  exceedingly  low. 
In  addition,  testing  electrically  is  much  less  time-consuming  as  the  steady  state  can  be 
reached  rapidly.  In  order  to  make  this  measurement  useful,  the  selection  of  test  materials 
is  critical.  The  specimens  should  not  be  chemically  active  at  room  temperature.  In  this 
regard,  good  candidates  include  gold  and  platinum. 
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MATHEMATICAL  EQUATIONS 


Probability  density  function, 

fix)  - 

Jin  a 


Error  function, 


Integration  of  relevant  exponential  function, 


/ 


,-Axe-Bx2x2(jx 


2 B 


xeAxeBx 


2 ib 


A 
2 B 


e-Axe-Bx2  |* 


l \ 

A2 

A2 

- - 1 

sfre  4B 

erf 

\ Zb  J 

2 s[B 

\[Bx  + 


A 

b 

2 s[B  ) 

a 

(A . 1 ) 


(A. 2) 


(A. 3) 


130 


APPENDIX  B 

RECORDED  SURFACE  PROFILES 
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a)  0° 


Vertical  scale:  2/xm/DIV  Horizontal  scale:  O.lcm/DIV 


b)  90° 


Vertical  scale:  2/xm/DIV  Horizontal  scale:  O.lcm/DIV 


Fig.  B.l 


Recorded  surface  profiles  of  Specimen  A:  a)  0°  and  b)  90° 
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Fig.  B.2 


Recorded  surface  profiles  of  Specimen  B:  a)  0°  and  b)  90° 
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a)  0° 


b)  90° 


Vertical  scale:  1/xm/DIV  Horizontal  scale:  O.lcm/DIV 


Vertical  scale:  1/xm/DIV  Horizontal  scale:  O.lcm/DIV 


Fig.  B.3 


Recorded  surface  profiles  of  Specimen  C:  a)  0°  and  b)  90° 
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a)  Oc 


Vertical  scale:  1/xm/DIV  Horizontal  scale:  O.lcm/DIV 


b)  90° 


Vertical  scale:  l/^m/DIV  Horizontal  scale:  0.  lcm/DIV 


Fig.  B.4 


Recorded  surface  profiles  of  Specimen  D:  a)  0°  and  b)  90° 
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